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liquid  b  I  props  I  lari  t  rocket  ure  derived.  Those  are  utilized  in  an  analytical 
Investigation  of  two  mechanisms  which  are  capable  of  producing  linear 
combustion  Instability  In  the  high  and  intermediate  frequency  ranges, 

typified  by  the  appearance  of  transverse  waves  and  entropy  waves 

1;,.  ‘f  '  ,V,  ■-  Hr  -  ..-.  x  J: 

respectively. 

The  character -1st ic  equation  of  each  rocket  system  is  derived, 
and  It  Is  shown  how  the  stability  limits  may  be  determined  for  a  particular 
rocket  motor. 
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1  *  introduction  to:  the,  theory. 

ir  -  / ■  't'fe  if/-/- r-'-rP  •' $  /  J ■  f  " ., 

k  iixlmteUim. 

Coupled  ’-j* 

over  I,  I SN-  ,or f dov«^SKf  (l ,  qh  |gf  g  ' 

chamfers,  there  lias 'dr  I  son  . ^  which  has 

SJrow  i»  ImporTonco.  Tho  nrohlnm  uki  i,  , 

p-ooiun,  which  car,  arjse  ,„  any  ducttHl  b|)r  , 

tho  phenomenon  of  combust |on  Instability,  and 

«nu  ■  l  llioucjh  a  proc  I  go  un  I  vorru.i  I  |  y 

-coptod  daMnltlon  does  not  exlsfy  It  Is  genera,  ,y  considered  o, 

r°9Ul<"'  P°ri0‘,IC  0SCM  9tl“*  «■*•»■*  chamber  press,.,,  o  wh  ,ch  arp  , 

fi’ln0d  ,M  mnmr  ,y  th9  — tton  process.  These  prWe  wrlaMw, 

flr°  0bSOr™d  U,<POrliBntdlly  10  Covor  ^  froquoncy  range  from  ,0  to  12,000 

CyC'9S  P“  SOCOnd' Whlcl1  WY  fr«  ^.ohe  hundi'od  percent 
of  cluimbor  prossurQ; 

°SCl ' l0tl0nS  °'  pre5SUrc”  a"d  S**  °[  «l  I  dependent  thermodynamic 
°'Kl  ,IUld  ClVn“"IC  Varlatl°3'  ara  '’lghly  u,l<los Irab |e  j„  general,  since  median I- 
°r  tll0r™al  ftlllUrS’  8nd  eon+ro1  malfunction  may  ensue  shortly.  The 
Wry"H'  Ch“nber  Pr°S5Ure’  '9S  M"  as  “hying  thrust  J,  unduly  stros,  or 
'T"9U°  ,i0th  th#  °tamb8r  °"[l  — •*  calls Ing  mechanical  failure.  Tho  ho#1 

t-nsfor  to  the  wall  can  ho  Increased  several  hundred  procent,  duo  to  the 

I  requeue  y  oscillations,  which  may  ho  sufficient  to  cause  rapid 
deter  loratlon  and  burnout  of  tho  chamber  wall.  P, no „y,  even  If  the  rochet 

“t0r  WithStand  th9  V,bra+I°"  *"«  ti ansf or,  severe  secondary  oscIMa- 

,,0"S  "y  1,0  *0t  “p  ln  HH  dolkaie  SUKtanco  and  control  system  which  will 

J"  '">y  ltS  0,tfBC'l'IVaneSS'  Af'V-  Oh  nil.  Of  these  effects  will  result  lu 

,a"Ure  Pr°PUlSi0n'  Unll',£,nd  0,n  |  from  consideration  for 

Inclusion  In  a  system  which  requires  rollahlllfy,  accuracy  arl(J  d()pondah , , ,  |y_ 

SI"C°  1110  Pr8r,0n°e  COmbUS,^£  effects  tho  life  and  re,|ah,|||y 

+h°  r°Cke+  m°t0r'  "  19  ,1C:'CGSGary  that  we  gain  an  understand,,,,,  of 


wh Ich  wl  1 1 


fluid  dynamic  processes  so  that  tho  conditions  which  promote  this  detrimental 
form  of  combustion  may  bo  dotermlnod  arid  consequently  avoided,  control  led- 'or 


■ft. 


•  A. 


el Imhia+od. 

Although  tho.ro  have  boon  other  attempts  at  classification,  (for 
oxamplo,  see  Refs.  15  and  II)  tho  frequency  spectrum  In  which  combustion  In¬ 
stability  occurs,  can  bo  con von  lent ly  separated  Into  three  parts,  denoted 
low,  Intermediate  and  high  frequency  Instab  1 1 Ity,  respectively.  This  op- 
pears  to  be  a  natural  division  since  the  coupled  wave  process,  governing 

■if 

ea>,h  of  theso  three  types,  are  different  and  hence  load  to  different  charac¬ 
teristic  frequencies. 


■If  ; 

1h: 


it 


If  all  the  gas  in  the  chamber  surges  periodically,  large  Inertias 
are  Involved  and  low  frequency  oscillations  result,  those  low  frequency  In¬ 
stabilities,  commonly  referred  to  as  "chugging",  have  frequencies  ranging 
from  10  to  200  cycles  per  second,  and  have  been  demonstrated  both  theoreti¬ 
cally  (Refs.  2  and  3)  and  experimentally  (Rot.  10)  to  depend  primarily  upon 
the  coupling  with  the  propellant  feedlines.  Since  this  case  has  already  been 
thoroughly  explored,  and  since  the  Instability  may  be  eliminated  In  most  In¬ 
stances  by  Increasing  tho  Injector  Impedance,  that  Is,  by  Increasing  the  In¬ 
jector  pressure  drop,  we  need  not  consider  this  type  further. 


intermediate  frequency  instabilities  have  boon  observed  experimen¬ 
tally  (Refs.  9  and  16)  and  generally  occur  at  frequencies  of  several  hundred 
cycles  per  second.  It  has  been  postulated  ( soo  Section  2)  that  the  character- 
Istk  frequency  may  be  attributed  to  the  presence  of  entropy  wavos  In  the  chamber. 
Mils  typo  h,as  not  yet  received  a  comprehensive  analytical  treatment,  until  now. 

High  fi  equoncy  or  "screaming"  oscillations  are  generally  associated 
with  various  acoustic  modes  of  tho  chamber  and  occur  at  frequencies  between 
several  hundred  and  several  thousand  cycles  per  socond  depending  on  tho  mode, 
chamber  geometry,  and  exhaust  nozzle.  That  Is,  experimental  observation 


Indicates  'that  the  intense  shrill  sound  that  occasionally  Issues  from  4 
^//ockot  .  Is:  produced  by  Itlgh frequency  .pressure  osc-l  I  lotions  character  Izod 
by  frequencies  In  the  neighborhood  of  the  organ-pl po  rosoiiaiicos  ol  tho 
chamber »  Those  In  stub  1 1  j  lies,  wore  Initially  found  to  bo  of  tho  longitudi¬ 
nal  type,  but  with  moro  adequate  Instrumentation,  It  was  discovered  that 
transverse  modes  wore  also  present.  Tho  most  frequently  encountered  Irens- 
vorso  modes  are  forms  of  the  first  tangential,  or  "spinning-sloshing"  mode, 

In  which  the  pressure  waves  propagate  diametrically  or  tangential ly  across 
Iho  chamber.  Transverse  modes  of  Instability  In  rocket  chambers  have  not  as 
yot  rocelvod  comprehensive  analytical  treatment,  until  now. 

In  all  of  those  cases,  If  oscillations  are  to  bo  maintained,  there 
must  be  some  process  of  coordination  which  periodically  foods  sufficient 
onorgy  Into  the  oscillating  gras  system  to  sustain  tho  process.  Of  course, 

In  order  to  have  a  closed  cyclo,  the  gas  dynamic  system  must  Induce  the  com- 
bus  Mori  process  to  release  energy  at  the  proper  time  phase-wise  during  each 
cyclo;  thus,  any  theoretical  model  must  demonstrate  a  closed  loop. 

Iho  analytical  Investigation  presented  hero  consists  of  two  parts. 
The  first  deals  with  tho  treatment  of  high  frequency  transverse  mode  Insta¬ 
bility,  In  which  the  first  tangential  mode  Is  of  primary  Interest.  The  second 
deals  with  a  theory  of  combustion  Instability  In  the  Intermediate  frequency 
range  which  Is  characterized  by  tho  appearance  of  entropy  waves  In  the  chamber, 
Although  analytical  treatment  of  tho  problem  of  combustion  Insta- 
bllliy  In  liquid  propellant  rockets  Is  comparatively  recent,  an  extensive 
I  iterature  on  iho  subject  already  exists.  Undoubtedly,  tho  most  comprehensive 
treatment  of  the  subject,  to  date,  Is  duo  to  Crocco  and  Chong.  Because  of  tho 
thoroughness  of  the  analysis  and  discussion  appearing  In  tho I r  rocent  treatise, 
(Rof.  IB),  tho  reader  Is  advised  that  tho  latter  constitutes  tho  primary 
reference  for  this  thesis. 
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In  their  recent,  survey,  (Ref.  13),  Putnam  and  Dennis  point  out  that 
of  all  typos  of  combustion  oscillations  which  have  boon  observed,  the  organ- 
pipe  type,  In  which  tho  wavelength  of  the  oscillation  Is  related  to  the  dimen¬ 
sions  of  tho  chamber,  has  the  oldest  history,  It  seems  that  Higgins  produced 
an  acoustic  oscillation  or  "singing  flame"  in  1777  by  surrounding  a  diffusion 
flume  with  a  largo  duct  open  at  both  ends. 

Since  that  time,  thermo-acoustic  oscillations  have  boon  observed 
In  many  different  pieces  of  laboratory  and  Industrial  equipment,  for  example, 
detonation  tubes  and  gas-fired  units,  and  consequently,  many  Investigators 
have  attempted  to  analyse  the  mechanisms  Involved, 

One  of  the  earliest  analyses,  due  to  Rayleigh  (Ref,  I),  was  one  In 
which  he  advanced  a  criterion  without  formal  proof,  that  for  the  excitation 
of  thermally  driven  oscillations,  there  must  he  a  fluctuating  host  release 
within  the  medium,  such*  that  It  has  a  component  In  phase  with  tho  varying 
component  of  the  pressure  at  the  position  of  heat  release.  Note  that  this 
Implies  that  heat  release  at  a  pressure  node  cannot  contribute  to  thermo- 
acoustic  oscillations,  In  a  crude  sense,  Rayleigh's  criterion  may  bo  con¬ 
sidered  to  bo  nothing  more  than  tho  restatement  of  tho  conditions  for  a 
closed  thermodynamic  cycle.  That  Is,  If  heat  is  supplied  to  a  medium  at  high 
pressure  and  rejected  at  low  pressure,  then  piston  work  may  bo  obtained. 
Presumably,  this  net  work  Is  utlllzod  In  driving  the  pressure  waves. 

As  previously  stated,  many  contributors  have  analysed  the  fluid 
dynamics  of  various  combustion  driven  oscillations,  for  example,  see  Refs.  12, 
13,  14,  IS  and  20.  However,  since  these  and  other  published  works  do  not  bear 
directly  on  the  specific  problem  which  concerns  us,  namely,  combustion  In¬ 
stability  In  liquid  propellant  rockets,  wo  will  now  proceed  to  review  some  of 
tho  concepts  loading  to  the  theory  presented  In  the  following  sections.  Wo  • 


will  begin  with  the  evolution  of  Iho  I  lino  lag  concept  In  rocket  motors. 

First,  we  recognize  that  the  normal  process  of  combustion  In  a 
rocket  motor  is  of  a  highly  turbulent  nature,  and  hence  there  are  time-  and 
spaco-wlso  fluctuations  of  'Iho  pressure  and  dependent  thermodynamic  and 
fluid-dynamic  variables  throughout  Iho  chamber.  If  the  fluctuations  have 
small  amplitudes,  then  the  combustion  Is  considered  smooth,  and  when  they 
have  large  amplitudes,  It  Is  called  rough  burning.  However,  these  terms  are 
qualitative  at  best,  inasmuch  as  I  he  amplitude  of  the  oscillations  alone  Is 
an  insufficient  criterion  for  classifying  tho  stability  of  the  combustion. 

Random  fluctuations,  l.o.  those  which  do  not  have  a  characteristic 
frequency,  exhibit  the  attributes  of  lurbulonco  and  In  certain  cases  will  not 
bo  detrimental  to  the  practical  operation  of  tho  rocket.  That  Is,  tho  time 
averaged  exciting  forces  arc  negligible  and  the  combustion  Is  rough  but  stable. 
When,  however,  per Iodic  fluctuations  are  present  such  that  one  or  more  com¬ 
ponents  of  tho  frequency  are  predominant  and  grow  In  amplitude,  the  Integrated 
effects  will  be  non-zero  and  mechanical  or  thermal  failure  may  follow  In  short 
order.  Accordingly,  wo  note  that  tho  occurence  of  rough  and  detrimental  com¬ 
bustion,  given  the  name  unstable  combustion,  Is  characterised  by  oscillations 
with  well-defined  frequencies  whose  amplitude  Is  limited  only  by  the  damping 
of  tho  system.  This  distinction  between  rough  but  stable  combustion,  and  un¬ 
stable  combustion  was  advanced  by  Crocco  several  years  ago. 

In  order  that  unstable  combustion  exist  as  defined  above,  some 
coordinating  Influence,  which  Is  capable  of  amplifying  a  random  disturbance, 
must  be  present.  In  this  connection,  In  1941,  von  Kerman's  group  advanced 
Iho  concept  of  a  combustion  time  lug,  or  delayed  Instantaneous  combustion. 

This  was  defined  as  tho  time  between  tho  Injection  of  a  propellant  element  and 
Its  evolution  Into  combustion  products,  which  was  assumed  to  occur  Instantane¬ 
ously,  after  a  certain  delay. 


Gur,der  a,l(J  fr  ,dnt>  (Rof*  2)  Worf)  the  first  to  present  a  formal  I  zed 
analysis  applying  tills  concept  to  a  treatment  of  low  frequency  combustion  In¬ 
stability.  Ttioy  postulated  a  constant  value  of  the  time  lag  and  were  able  to 
demonstrate  that  under  certain  circumstances,  tills  time  lag  could  provide  the 
necessary  coupling  between  I  tie  rocket  chamber  and  the  propellant  food  system 
which  would  result  In  unstable  amplification  of  a  random  pressure  disturbance. 
Other,  more  elaborate,  low  froquoncy  analyses  followed,  Including  the  early 
work  of  Summer f told  (Ref.  3)  who  established  the  theoretical  limits  of  low 
froquoncy  stability  In  Ills  analysis  wlilcli  Included  the  effects  of  chamber 
capacitance  and  feeding  system  Inertia.  However,  since  each  of  those  retained 
the  constant  time  lag  concept,  the  variation  of  the  rate  of  burned  gas;  genera- 
"on,  or  source  of  driving  energy,  could  bo  modified  only  by  the  rate  of  supply 
of  propellant  to  the  chamber,  which  dopendod  entirely  upon  the  sensitivity  of 
the  feeding  system  to  chamber  pressure  oscillations. 

Now  although  Crocco  was  willing  to  retain  the  concept  of  an  Instan¬ 
taneous  transformation  of  liquid  propellant  Into  gaseous  products,  ho  reasoned 
lhat  the  specific  rate  of  conversion  of  unburned  propellant  elements  Into  com™ 
bUSt l°n  Procluc+s  depends  on  the  sensitivity  of  the  activation  processes  to 
chamber  oscillations,  and  In  |9I>|  he  Introduced  the  concept  of  the  sensitive 
time  lag.  Thus  the  total  time  lag  Tf  ,  during  which  various  physlco- 
chomlcal  processes  occur,  was  considered  to  consist  of  two  parts,  a  constant 
port  ^  durln9  whlch  mechanical  processes  Insensitive  to  the  thermodynamic 
states  of  the  surrounding  gas  take  place,  and  a  variable  sensitive  part  T 
dining  which  activation  processes  occur  which  are  sensitive  to  the  oscillations 
occurring  In  the  surrounding  gas,  Simply  wrltton: 

Tt  .  %  +  r 
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whoro  the  bar  denotes  s+oody  stole.  And  then  to  relate  the  sensitive  lime 
log  lo  the  rates  ol  the  conditioning  processes,  Crocco  wrolo  on  Integral 

oq unt  Ion  for  an  olomonl  burn  I  ncj  o  l  Nino 
/ 1* 

j,  “f  (i1)  d{ '  a  consf, 

't-r  2.2 

whoro  f  I u  i ho  ovoro 1 1  r.rlo  of  tho  conditioning  processes,  ond  tho  Integral 

musf  bo  evaluated  following  oach  Individual  propellant  olomonl.  Tho  rato  *F 
vorloo  along  tho  path  of  Integration,  and  upon  corrolating  all  of  tho  physi¬ 
cal  factors  which  cause  this  variation,  to  the  pressure,  Crocco  obtained  tho 
re  I  at Ion 

f s  f  ( 1  +  'ft-  ■§■ ) 
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whli.li  Is  valid  for  small  amplitudes  of  oscillation.  Tho  quantity  \  which 
Is  actually  a  well-defined  mathematical  quantity  (soo  Ref.  10)  was  called  an 
Interaction  Index,  and  was  assumed  to  ho  a  characteristic  of  a  given  propellant 
combination.  Crocco  observed  that  actually,  of  course,  sensitive  and  Insensi¬ 
tive  time  lags  are  physically  Inseparable  and  occur  simultaneously.  Hence 
Tq.  ?..|  should  be  considered  os  a  schematic  representation  for  the  actual 
events.  On  combining  Eqs.  2.2  and  2.3  we  observe  that  tho  sensitive  ! Imo  lag 
exhibits  Inverse  pressure  dependence,  since  there  Is  readily  obtained: 


T  p  “  Con  si. 
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Application  of  tho  sons  It  I vo  time  log  concept  to  tho  treatment  of 
combustion  Instability  was  first  mode  by  Crocco  (Ref.  4)  for  tho  cases  of  low 
frequency  Instability  In  monopropol lanl  and  b I  propel  lent  motors,  and  roughly, 
for  longitudinal  high  froquoncy  Instability.  Tho  term  Intrinsic  Instability 
was  Introduced  to  Indicate  that  duo  to  tho  coupling  between  the  pressure  oscll- 
lal Ions  and  combustion  processes,  through  tho  medium  of  tho  sensitive  time  lag, 
Instability  can  occur  oven  if  the  Injection  system  delivers  a  constant  flow. 


a 


A  series  of  papers  l.y  Crocco  and  Cheng  dealing  with  n  more  refined 
treatment  of  the  low  and  high  frequency  longitudinal  modes  of  combust  Ion  In- 
sfablll'ly  followed  and  Was  recently  assembled  In  their  monograph  (Ref .  |H). 

Ill  Mils  same  reference,  It  was  pointed  out  that  theoretical  analyses  wore 
lacking  for  two  Important  mechanisms.  Accordingly,  wo  will  consider  the  treat 

"lent  of  those  two  distinct  problems,  transverse  wave  and  entropy  wave  com- 
bustlori  Instability, 

In  both  casus  to  bo  analyzod  horn,  the  troatmoirt  Is  restricted  to 
llnoar  Instability,  whore  tho  perturbations  are  sufficiently  small  so  that 
second  order  terms  may  bo  considered  negligible.  It  Is  recognized  that  fully- 
dovolopod  combust  lor.  Instability  Is  general  ly  a  nori-l  Inoar  phenomenon, 
characterised  by  the  pretence  of  shock  waves,  but  as  has  boon  observed  ox- 

|,or  h"°"'l'ally  ln  0  9roa+  manV  cases,  Instability  will  oftlmes  result  from  tho 
pi  ogress  I vq  amplification  of  small  disturbances.  Thus,  when  wt  determine 
stab!  "  I y  limits,  It  moans  that  within  those  limits,  tho  system  will  bo  stable 
to  small  disturbances.  If  a  giver,  system  Is  linearly  stable,  and  no  largo  dis¬ 
turbances  aro  applied,  then  self-ampllflcatlon  cannot  drive  tho  system  to  In¬ 
stability,  and  no  hlgh-ampl Itude  pressure  oscillations  can  appear. 

Lot  us  now  describe  tho  physical  system  and  tho  mochanlsms  which 
are  Involved.  At  flier  high  frequencies  normally  encountered  In  transverse  mode 
Instability,  the  Injection  system  cannot  respond  to  chamber  pressure  oscil¬ 
lations.  Thus,  for  Instability  to  exist,  It  must  be  of  the  Intrinsic  variety, 
and  a  coordinating  process  must  ho  „rosont,  so  that  oscillations  of  tho  rate 
I n g  factors  will  produce  organized  oscillations  of  the  burning  rates, 

which  will  provide  the  necessary  exciting  force  to  maintain  the  coordinating 
process. 

Iri  out  Investigation  of  transverse  modes,  wo  will  retain  Crocco' s 
m^dol  of  tho  sensitive  time  lag  as  a  suitable  mechanism,  and  wo  will 


Investigate  tho  stab  i  n  f  y  of  three-dimensional  perturbations  with  exponential 
tlifio  dependence,  In  a  general  rochet  r.yntom  cons  luting  of  n  rigid  cylindrical 
chamber  with  a  liquid  propel  I  an  I  Injector  and  a  fixed  convorglng-dl verging 
exhaust  nozzle.  It  should  bo  observed  that  duo  to  the  presence  of  tho  exhaust 
no/zlo,  purely  transverso  oscillations  cannot  bo  present.  If  a  transverse 
perturbation  appears,  a  longitudinal  wave  will  be  rof  looted  in  the  subsonic 
portion  of  111 o  nozzle,  and  thus  a  form  of  combined  mode  must  exist,  In  contra¬ 
distinction  to  what  occurs  in  a  cylindrical  chamber  terminated  by  plane 
closed  ends. 

.  To  explain  the  experimentally  observed  occurrence  of  combustion  In¬ 

stability  In  the  Intermediate  frequency  range  (Ret.  9),  another  mechanism  must 
be  considered.  In  Mils  case,  tho  frequencies  are  too  high  to  ho  the  result  of 
chamber  and  feed  system  coupling,  and  they  are  too  low  to  correspond  to  a 
mod©  of  resonance  of  tho  gases  In  tho  chamber*  Thu  proposed  mechanism  depends 
on  the  presence  of  entropy  waves  In  I  ho  chamber,  and  In  some  cases  involves  a 
special  form  of  coupling  between  the  Injection  system  and  the  chamber. 

Entropy  waves  may  be  produced  by  chamber  pressure  oscillations  in 
one  of  two  ways,  each  of  which  Is  capable  of  reinforcing  tho  other.  Howovor, 

In  a  b I  propot  Ian  f  rocket,  tho  primary  cause  arises  when  one  obtains  a  varia¬ 
tion  of  tho  mixture  ratio  about  Its  moan  value.  This  occurs  when  Iho  oxidizer 
and  fuel  injectors  respond  differently  to  chamber  pressure  oscillations,  thus 
producing  a  stream  of  propellant  elements  with  an  off -rn Ixturo  ratio.  If  tho 
mixture  ratio  Is  oscillating,  tho  temperature  and  tho  entropy  of  the  products 
of  combustion  at  a  given  station  will  also  oscillate.  For  typical  rocket’s 
designed  tor  maximum  thrust,  both  become  ,u.  w,  or  smaller  as  tho  oxldlzor- 
fuol  ratio  Increases  or  decreases.  (sooFIg.  2.1) 

As  a  result,  If  one  examines  the  conditions  In  the  chamber  of  Ibis 
rocket  at  any  given  Instant  of  time,  one  observes  that  +>  o  gas  entropy 


10. 


cl  loir  I  but  Ion  has  a  wavo-liko  pal  torn  which  moves  downstream  with  ‘the  gas  veto- 
city.  The  entropy  oscl I  let  Ion  re I  loots  pr  assure  waves  at  I  ho  exhaust  nozzle 
which  travel  upstream  with  the  speed  of  sound  to  the  Injector  face,.  This 
closes  the  feedback  loop  allowing  the  process  to  start  all  over  again#  i he 
explanation  of  this  mechanism  was  first  given  roughly  by  Dorman  arid  Cheney 
(Rof.  9),  and  has  since  been  enlarged  upon  (Refs.  16,  10).,  In  the  mathemati¬ 
cal  Iroatinont  In  which  the  stability  limits  for  this  mechanism  will  bo  deter¬ 
mined,  Crocco's  model  of  the  sensitive  time  lag  Is  ext  ended  to  allow  for  the 
functional  dependence  of  the  rate  of  the  conditioning  processes  on  the  mixture 
ratio.  The  treatment  of  entropy  wave  Instability  appears  in  Sections  9 
through  1 3. 


II.  TRANS  VERS!  WAVE  INSTABILITY 
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In  a  goriora 1 1 zed  iruatmenl  of  combustion  lnstab 1 1  Ity,  wo  are  con¬ 
cerned  with  I  ho  throo-d  Imonr.  Iona  I  motion  of  a  cjao,  containing  a  dish  Ibutlon 
of  droplets  of  liquid  propellant,  which  flows  Through  a  rockot  chamber,  Tho 
droplets  which  burn  at  different  locutions  In  tho  chamber  will  then  corres¬ 
pond  to  0  distribution  of  sources  of  mass,,  momontum  and  energy.  Lot  us 
formulate  tho  conservation  laws  for  this  two- phase  flow. 

Wo  lot  P  denote  tho  density  of  the  gas,  defined  as  tho  mass  of 
gas  per  unit  volume  of  gas,  and  lot  ^  denote  tho  density  of  tho  liquid 
drop  lots,  defined  as  tho  mass  of  liquid  per  unit  volume  of  gas,  that  Is,  we 

,■  1  .11  ■:{ 

will  neglect  tho  volume  of  tho  drop  I  ots  with  rospocf  to  the  gaseous  volume. 
Then  lotting  a  donoto  a  dimensional  quantity,  tho  conservation  of  mass 
bocomon 


ip  ffr7f)  +  v*if^+0)  =° 
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where  \l  *  and  V*  are  the  gas  velocity  and  liquid  drop  I  el  velocity  ro- 
-6" 


I  rvk 

.poc lively,  't  Is  the  time  and  V-  Is  tho  dlvorgenco  operator, 


way 


Now  note  that  tho  continuity  equation  may  bo  treated  In  another 

■  t  1 

<q>  denotes  the  Instantaneous  rate  per  mill  volume  a  I  wh  h1'  gas  Is 


I  ir 

generated  al  any  local  Ion  In  tho  chamber,  so  that  Sr  corresponds  lo  a 
source  distribution  for  the  gas  phase  and  simultaneously  a  sink  for  tho 

liquid  phase,  then  continuity  may  also  ho  written: 

4 
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If  viscosity  Is  neglected,  itio  conservation  ol  momentum  In  this 


two-phoso  system  takes  I  ho  form: 


which  Is  sornowhat  more  comp  I  I  Ctrl  oil  Ilian  I  ho  faml  I  lor  In  lor  equation. 

Tho  energy  changos,  must  obey  the  first  law  of  thermodynamics, 

I . o . ,  tho  wor  k  done  on  tho  gas,  plus  tho  boat  added  to  tho  gar.  must  equal 
tho  change  of  energy  of  tho  gas.  The  several  forms  of  energy  which  will  bo 
considered  bore  Include  kinetic,  Internal  and  chemical  energy,  however,  I  ho 
work  done  by  viscous  stresses  and  heat  transferred  by  conduction  or  diffusion 
will  ho  neglected.  13  y  dof  In  1 1' ion, 


f 
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and  noting  that 


i  ^ 


since  for  tho  liquid  phase  I  lie  Internal  onergy  and  tho  enthalpy  are  very 

l  ^ 

nearly  tho  same,  and  tho  corresponding  common  value  V\j  Is  Intended  to  In¬ 
clude  the  chemical  energy  of  tho  propel  Inn  It.,  we  obtain  the  following  equiva¬ 
lent  forms  of  the  energy  equation: 


13. 
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Not  I ng  that  n  Isa  function  of  the  temperature,  the  four  equations 
3.1,  3.2,  3.3,  and  3.7  contain  eight  unknowns  ,  N/  j  p  ■>  I  j  T  ■>  ,  VjL 

and  \)£  ,  lionco  four  additional  equations  are  required.  First  wo  have  the 


equation  of  state  for  t ho  gas  phasos 

4  ..  />  *  ro  VI  T  + 
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A  second  equation  relating  the  burning  rate 
quant  1 1 1  os  will  bo  dor i vod  I  at  or . 


to  the  other 


A  third  equal  Ion  Is  obtained  from  the  dynamic  behavior  ol  the 

drop  I ots : 

.  ’I*  Ml*-V) 
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which  assumes  that  tho  forco  oxortod  by  the  gases  on  the  liquid  droplets  Is 

Inversely  proportional  to  tho  Reynolds  Number. 

Tho  fourth  equation  could  he  obtained  from  the  heat  balance 

I  4  .  f- 

of  tho  droplets  which  would  yield  mi  expression  relating  v\^  and  n 
We  shall,  however,  neglect  1  ho  boat  transfer  ho  I  ween  the  gasos  and  I  ho  liquid 


droplets  by  faking  (sool'tof.  IB) 
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This  Implies  Hurl  when  wo  (allow  1  ho  mol  Ion  sf  a  pai  I  icular  droplet,  the 
\  ^ 

value  of  Is  conserved.  In  other  words,  tho  droplet  iota  I  nr.  tho  value 

\  # 

of  m  t  w I  I h  which  It  was  Injected.  We  may  then  wr  ite: 


\  #■  \  4  \lA%  \  %  .  4"  v/  ff2- 

Hi  ■-  "j  +  i  -  ty,,  a  In^,  +  Ve  3.11 

2»  2. 

\  \hl  ^  ^ 

where  niji0  and  are  t|io  values  of  the  propellant  enthalpy  and  kinetic 


energy  at  the  Injector  head,  and  we  remark  that  those  quant I  ties  may  bo  ro™ 
la  hod  to  tho  pressure  at  tho  Injector  by  means  of  tho  injector  response  aqua' 


I  Ions. 


Iho  Injection  system  does  not  respond  lo  chamber  pressure 


M. 

■Ju 

oscillations,  1  he  itiixl  m  e  i  alio  remains  eons  I  an  I,  therefore  111 L1 s  I  I  I  k  o  - 

\jf  *  *- 

wise  bo  ooir.il.ml,  and  unions  I  ho  inject  Ion  velocity  lb  modulated,  In 

.also  constant.  In  gr  nor  a  I,  llu-  rnlxturo  ratio  and  injocl  Ion  volot  Ity  are  not 
coir. tan  I,  and  each  droplet  rot  a  I  ns  n  different  characteristic  stagnation  en¬ 
thalpy  right  'through  the  moment  of  Its  conversion  into  burned  gas. 


Wo  may  now  prococxl  to  non-dimensional  I zo  llio  tor  ego  lug  equations 

**"'  *****  ^,i  - - —  1  »■  11  " 

by  using  ns  reference  values  ?S,  f\*  ,  To  arid  to  ~  ^  yCR^To^  , 

which  are  the  steady  stale  values  of  p^csu  re,  dons  II  y,  tempera turo  ol  the 

gas  and  the  voloc  Ity  of  sound  In  the  gas  at  Ilia  Injector  face.  Note  I  ha  I 

7)  * 

since  V_  must  be  zero  at  the  injector  faco,  the  rot  or once  values  corres¬ 
pond  lo  stagnation  values  of  Iho  respect Ivo  quant  It  Ion.  In  I  he  study  ol 
transverse  waves,  it  Is  convenient  lo  use  1  lie  chamber  radius  Ic  as  the 
reference  length,  whereas  In  Iho  study  of  longitudinal  waves,  the  comb us  I  ion 
chamber  length  Is  a  more  suitable  reference  length.  Wo  will  begin  by  select¬ 


ing  ft  as  the  measure  of  length. 


whore  y  is  the  ratio  of  specific  heats  which  Is  assumed  to  bo  constant  over 
the  range  of  variation  of  Hf  *  *  This  assumption  yields  din  a  d  i  I n  I  or ms 
of  the  non-dimensional  vai  lab  I  os.  The  non  dimensional  equations  may  now  bo 
wr 1 1 1 on : 


X  r)n 


1 

51: 


Y)l(v-v)  +  (Y-vyv  4 ^i(v.h) 
4  (Vt-v)^  _-  -i  vp 


3,3a 


( Y-v)kj 


ti  22- 
* 


4>(hj“hfj)  3.7ti 


P  B  z5 


+•  (Vt-vjVe 


3,  Ua 


3»9a 


ss  O 


V^S  =  S 


,3 1, 1  On 


I  Id 


and  now  lot  us  pi  oeoocl  by  introducing  small  porturbat Ions. 

Wo  will  cons  I  dor  each  of  I  ho  dopmidont  variables  as  I  ho  sum  of  a 
steady  state  spaco-var lab lo  and  a  I Ime-dopondont  porturbat Ion  so  small  that 
forms  higher  than  those  linear  In  the  perturbations  can  bo  neglected.  Thus 
p  =  p  f  p^  ?  p  •-  p -H  p1  ole.,  whore  the  superposed  bar  denotes  steady 
state  and  the  prime  denotes  small  per  fur bu I  Ion.  Tho  steady  slate  equations 
and  I  he  equations  linear  In  tho  porturhat Ions  follow  directly. 


16. 


Conti  mi  j I y 


17, 


(v.  7)% 


S  v(%.v)v'  +  (v.v)yi  ^(V-V) 


Heat  Trans for 


=  K 

y  -  V) 


/ 


r-\ 


Y  V  *  Ku 

(y"0  Ve, ‘V^1  s  in 


jj2$o 


5,9b 


3,9c 


3. 1 1  b 


3,  I  Ic: 


Al  Hi  In  point,  we  have  obtained  a  sot  of  pnrl  lal  differential 
Ch|ua'l'lon=  flovorn  the  motion  of  the  rock.. I  fluid  system  during  steady 
and  unsteady  operation.  Wo  are  now  Interested  In  determining  the  stability 
of  solutions  with  exponent  In  I  time  dependence,  hut  first  wo  proceed  with  the 
selection  of  an  appropriate)  coordinate  system. 

4* 


I  ho  nufufeil  coordinate  system  for  rocket  motors  of  conventional 
shape  Is  a  cylindrical  coordinate  system  and  hence  tho  latter  will  be  utl- 

llMd  horu-  Roforrln9  t0  rig.  4.1,  wo  see  that  I  no  square  of  the  elementary 
length  Is  given  by: 


a  1*' 

Non  dimensional  U I  rig, 
Z  * 

we  obtain 


de*1"  +•  dr**'2  + 


y-  63. 


rc* 


a 


4.3 


which  enables  us  lo  expand  I  ho  vector  equal  Ions  we  have  derived  previously. 
In  the  ensuing  inn  lysis,  It  will  ho  assumed  Ilia  I  the  steady  state  solution 
consist's  of  ono-d  I  mens  Iona  I  (low,  and  then  the  so  I  ill  ion  In  unsteady  flow 
will  consist  of  three  dimensional  perturbations  superposed  on  the  steady 
stale  solution.  It  follows  that  all  partial  derivatives  of  steady  state 
quantities  with  respect  lo  f  or  0  vanish,  while  partial  dorlvatlvos  ol 
steady  state  quantities  with  respect  to  g  become  ordinary  derivatives. 

Lotting  the  subscripts  £• ,  Y'  and  0  denote  components  ol 
vectors  In  the  respective  directions,  the  continuity  equations  become: 


di  (  fa' 


Proceeding,  we  obtain  the  equations  for  the  conservation  of 

momentum: 


19, 


"uwr  '"0  i no  onorgy  equation, 


V4M  vv»  ! /I  I 


—  IIIMVJM  \  \  UUO  arifl  I  yb  I 

5h0WS  th0t  h»  18  Constflllt  throughout  the  chamber  to  within  forms  of 

order  of  tho  squoro  of  the,  local  Mach  number  and  honco  from  Eq.  3.7b  If 

111,11  ^  *  Vt  '  Wo  1,111  J^Ufy  this  ordor  of  mngnltudo  anala- 

11  lG  ,n  another  soctlon.  Tlioroforo, 

C\  j-  n  V n  r)h.v  v  i 


*  t*  $  -  £  (w~bjs) 

»iay  $  by  Introducing  Eq.  4,4  and  we  can  ^Imlnafa  Kj  as 

fol  lows,  since  dh  =dT  ,  WO  must  also  have  h‘  =  T  '  and  therefore 

hs  =  T'  +  (ir-i)  V-V1  =  t1  .  ,  .  w,  u' >. 

-&>  I  •  ( ^ • )  V^S  Vz»  •  „  Introducing 

llio  oquatlon  of  state,  T  1  ~  T  /  P1  /) *  \ 

\  p  W  )  Wo  0,J,7llri  for  tho  onorgy 
oquatlon:  '  r 

a|  P'-TP1  +(jr-i)yO  V4 Vi1  1+  2.  f  c/a  .\->7  „ii  9 


j^^P'-T/o'  +  M^W}  j 


o' -I 

y  r 

0  s-t 


+  'Vs  ii 
1  da 


20. 


Iho  dr  op  1  o  l  dyn<: 

miles  yield: 

DjU 

V  i.e-  r 

dl  ( 

A  ( -v«t ) 

4.11 

+  %  1^4 

Si 

+• 

J  ii 

~  A  (4  ) 

4,12 

2Vr 

H 

4  s\^r 

s* 

=  Jk  (Vr1 

-Vt!-) 

4.13 

+  Ve*  aft* 

3i 

=  A  ( Vi? 1 . 

™  ^4$ ) 

4.14 

whllo  tho  droplet  stagnation  enthalpy  Is  b>p|y 

^  +  (tVjUt  s  h^u 

•j*  SPJMQlt  Ion  Of  the  for  jail  I  rn; 


Wo  may  aoparafo  tho  variables  In  our  system  of  linear  partis  dif¬ 
ferential  equations  by  exproaslng  each  of  tho  perturbations  aa  tho  product 

°f  . . tlono  0#  lhH  coordinates  and  an  exponential  time  function.  After  «x~ 

amlnlng  purely  arbitrary  functions,  It  Is  found  that  separation  may  be  achiov. 


who n  wo  take: 


v.' 

Vr 1 
Vo' 


V. 


% 


(e)  V(r)  $  (&)es^ 


■pr  (i)  MO  $  (d)e;i 


5. 


5-fc 


W  =  *)«  (*)  W  $  (o)  ed 

Vir1  =•  y  (a)  5  (o)est 

Yto' 53  ")o  (»)  iM  i0L(e)eii 

a'  =  J(6jest 

p‘  =  P(&)  ft*-)  §(0)est 

t'  =  l(1')  5  (e)e^ 
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S I  nee  1  ho  exponent  $  *  JV.+ L$ l  I  s  genera  1 1  y  comp  I  ox,  I  he  stab  I  f  j  ty  o  f  any 
solution  will  depend  on  A.  .  That  Is,  the  solution  will  bo  stable,,  neutral 
or  unstable  as  Jt  Is  loss  than,  equal  to,  or  exceeds,  zero,  We  will  deter¬ 
mine  the  conditions  yielding  neutral  stability,  Inasmuch  as  a  knowledge  of  The 
stability  boundary  (  A.  «  0)  and  the  unstable  side  will  bo  sufficient  for 


our  purposes. 

Wo  note  that  although  the  perturbations  are  given  as  comp  I  ox  quan¬ 
tities  for  the  purpose  of  Indicating  phase  rotations,  only  the  real  parts  have 
physical  moaning.  Before  solving  the  perturbed  equations  wo  refer  to  Ref.  ID 
for  the  solutions  to  the  steady  state  equations.  Steady  state  continuity 
and  momentum  yield  . 


f 


V 


^  (? )  <r 


5.2 


rr 


,0  ,  ,  .  =■  wr  =  ficV let  -  fi% 

whore  represents  the  known  propollarit  injoctlon  rate  and  uT  is 

s 

the  gas  flow  rate  at  any  station,  and 

p  =  yA’ 

wli  1 1  a  wo  may  also  write 

7  «  I  -*■ 


1 

Utilizing  Eq.  3.0b  wo  also  obtain: 

_E 


V. 


a 


3.3 


5.4 


r 


3.5 


I 


% 


i-(y-i)  V* 

Considering  orders  of  magnitude,  wo  note  that  the  velocity 
Increases  from  zoro  at  the  Injector  aid  to  a  max’mum  value  connected  with 
The  Mach  number  at  1  ho  entrance  to  the  nozzle.  A  solution  of  the  non-linear 
differential  Equation  4.11  (see  Mg.  4,2)  shows  that  droplet  velocity 
'Vi  has  the  same  order  of  magnitude  as  V&  and  .  That  Is,  V*£  , 
and  ft  are  each  of  the  order  of  the  Mach  number,  furthermore,  the 
deviations  of  p  ,  ^  and  T  from  unity  are  of  0(M  ),  Up  to  twi  ms  of  this 


22. 


order,  wo  can  therefore  wrlto: 

p  e  f  =  T  =1 

and  honco  there  follows  from  Cq.  5.2. 

PM  - 

sinco  combustion  Is  assumed  complete  at  tho  exit  of  the  chamber  and 

m-rr- 

hi  Vti  ■  ■ 

If  wo  now  toko 

I 0  Tr  tSa  oft) 


If  wo  now  take 

s  o(0 

1 

% 

i&f 

d-fc 

Jfc 

\ 

Svr 

S  o(M) 

% 

1 

% 

^  §o(0 

and  If  the  Injector  does  not  respond  to  chamber  pressure  fluctuations 
(as  In  a  study  of  Intrinsic  Instability),  then  the  droplet  Equations  4.12, 
4. 13  and  4, 14  yield  simply: 

M  2.  ‘  d  I  l/\  m  ^4  I  ,A  So  i  .  5.9 


V  f  ^  f  s  f  ^ 

It  will  ho  shown,  In  our  treatment  of  entropy  wave  Instability  that  a  more 
complicated  result  Is  obtained  at  lower  frequencies,  since  tho  Injection 
system  can  then  respond  to  chamber  pressure  oscillations. 

And  now  substituting  from  Eqs.  5.1  Into  tho  continuity  Equation 
4.5,  wo  find  for  the  loft  hand  side: 

- — - - y  ><r  r’-t 


sS  +  +  i  (^S)“ 

PdJr 


d’-s 

dG*-  5J0 


whore  the  var lab  I  os  are  not  separated  as  yet.  We  will  return  to  this 
equation  presently. 

Upon  In+roduc  Eqs.  5.1  and  5.9  Into  the  equation  for  the 


23. 


component  of  momentum;  t-q<  4.7;  and  rearranging.''.  torms,  m  ob lain: 
cl 


+ 


fS 


5. 


i  if 
If  ii 


>  +  4S  V"'/3^r 


7  V? 


+  l/4$  +  'wi*\ 


d'Y  <W  „  J 

I5_  +  JE  +  ifc  Jli 

^  rf  vr 


Thu  T  and  d  components  of  momentum,  Eqs.  4.0  and  4.9  yield: 


A 

(1?: 


(pv*  + 1 


4*  S 


khM» 

f 


-I-  }if 


V3 


3.  12 


and 


+  S 


f  +  J/J& 


M 


9 


1 

if 


5.13 


Expanding  Eqs.  5.12  and  5*13  and  then  subtracting  one?  from  tho  other  m 
obtain: 


Lf  Vt+S 


Tills  has  tho  solution: 


ty-VAaO  s.K 


=  C,  exp  (-  ( 1  + 


5»  I  5 

i  -  s'('i  *Ji&  j  ’■  “  ■  "  v 

f  Vi  -I-  -fe  o. 

and  hone©  unless  G|  is  tnkon  Idontlcal  ly  ;:oro>  separation  of  tho  variables 
Is  procludod. 

II  Is  of  some  Interest  to  soo  what  happens  to  tho  vortlclty  as  a 
consequence  of  this  last  result.  By  definition;  tho  vortlclty  Is  given  by: 


rotV*  =  yxv'  -  4 

—+*  ■  -$»  i 


r 


•p 

***** 

30 


(  ]  Ui 

J  Hf 


5.16 


>(*’'  -  Ur  -  &M]  a 


i  m 


"  H 

i  . 


~-i  ...  -*r 
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Substituting  from  Cq,  5.1  wo  obtain 

(r0+l')a  =  ^  ~eSt 

(ro+v')r  ”  pit 

s  '  J  n.  /  -*«—* «  C_ 

r  ae 

•  (&-*,)£  i*rt 


2  O 


and  thus  wo  soo  that  the  condition  that  enables  us  to  separate  tho  variables, 
name  I y  Pg  'Vr  ,  also  causes  tho  axial  component  of  tho  vort Icily  to  vanish. 
Returning  to  continuity,  Gq.  5.10,  and  Introducing  p$*pr  >  wo 


f  Ind: 


pPr 


|rl.j'  yh  (jj 


wtier®  S  Is  the  separation  constant,  and  so  also  obtain 


df 

l  W.MWIIMWI 

+  1  dr 


H~  <5  nh  * 


whore  D  Is  the  second  separation  constant.  We  have  thus  obtained  the  set: 
f  ^  Jr  +  nh  ~  VV)  -O  5,20 


f  fv 


Equation  5.20  Is  recognized  as  a  Bos so  I  equation  so  that  I  ho 

solutions  become  respect  I vo I y, 

m  _  5  Jn(s«Var) 

1  Vva(Snhr) 

5  a  /  C.0S  n(? 

)  sirt  n9 


No  constants  nood  appeal  above  as  I  hoy  may  by  conveniently  ab¬ 
sorbed  Into  tho  *-  dependent  factors.  I  ho  actual  solution  Is  obtained  by 
summing  over  all  combinations  of  products  of  solutions.  Wo  must  satisfy 
corta In  conditions  Imposed  on  Vr'  .  First,  Vy'  must  not  become  Infinite 


at  ^  °>  0luJ  lu,ncH  «"  "'“St  discard  Hio  second  Bor, sol  function.  Further, 

Vr  mu  si  vanish  ah  Iho  cylinder  wall  where  f  |,  and  honco  Eq.  5. | 


Imp  1 1  os  that: 

dt\ 

I 

df  /r » 1 

-  dT  u«v)  *  o 

dr 

5.24 

Slnco  tlioro  aro  a  doubly  Infinite  number  of  solutions  to  Ibis  equation,  wo 
let  SX  represent  tho  h1h  zoro  of  tho  derivative  ol  tho  Bossol  function  of 

order  n,  so  that  both  1 

subscripts  of  SvnVo  aro  used 

as  Indices  to  Indicate 

which  of  tho  solutions 

Interest  us.  Tho  first  nine 

zeros  aro: 

S o\  m  0 

$02.  m  3,8317 

So^  -  7*0156  3*25 

Sit  *  1.8413 

Sri  *  5,3313 

S13  m  8,5263 

$a\  163  3,0543 

Su  m  6,7060 

$23  «  9,9695 

Tils  system  of  equations  which  must  ho  solved  simultaneously  now 

'ocomos ; 


&D±lmi±y. 


v 


sS  +  1  (p  44  +  V*  & )  -  = 

I  (  +  %  S)  + 

amato  * . 

+  Va£)  +■  +  Vj1  S)  =■  ^  ( 

f v*  ~ 5  ( 'if  Pa  +  Via  t ) 


5.26 


Snh 

d 


J  %  i 


T*  (°’  f  fa  ^  3 ) 


5.27 


26. 


Momoi  1 1  um  r 

s(^V  +  j-(f 


Easm 


v?  |  r -  T 


S 


</> 

If 


r 


T X  +  (y~i)^  V*24  no  ;'29 


Lot  us  rovlow  tho  Ini  Mol  conditions  at  tho  Injoctor  face  for  I  ho 


caso.of  Intrinsic  Instability  studios.  At  »  0, 

V*'  (  0, r, 0,-t)  =0  5i 

/V  (,0,r  0(-t)  s  0 

4>'  J  =  O 

Vfc  B  ,i  )  s  O 

slnco  at  tho  Injoctor  face  thoro  Is  no  production  of  go  3,  and  sjneo  tho 
liquid  is  Incompressible  thoro  can  bo  no  variation  In  propoll  an  I  density. 


Equivalently,  wo  have 

P%  (o)  ~  ^  (o')  s  J  (o)  a  ^Jp(o)  "  Vfc(o)*0  Jb”j| 

while  ^(0)  .  yg  ,  P^(p)epro  ,  and  £(0)  3  £o  . 

Wo  will  proceed  to  solve  oui  system  of  equations  by  an  Iteration 


scheme,  (soo  Ref.  18).  Combining  contlnu’ty  and  energy,  Eqs,  '3.26  and  5.29, 


and  lotting 


X  (%)  *  -t*  O-t) 

V  (i )  * 


wo  obtain: 


^  ^  4  "  ^  V  ^  5-  32 


d* 


§ )  +  i /  -J  \  -sk  =  -sX  +  £t  -sVk V 

‘ft  f  l  V  ^  /  ‘ft  ell 


5.33 


Reaming  I  rtg 

d  /  'f  \ 

d'i  {  W.  , 


1*  !)«-  /  w 0  ob  I  n  I  n : 

^  -  p' 

ilUi  < 


i-  s  t!i. 

v?. 


-{•  s^h  u 
j* 


5,34 


where  wo  have  token: 


wt? "  *&%+**$ 

u/o  "  3V  V  +  v~i  *v,;» 

u(*)  0  />V*70  +  *  W  * 

Out  last  equation,  I  q.  3.20  becomes:  ™ 


3.33 


s35r  +  Jt 


S  |Vl  -  olu 


5.36 


upon  taking 


M(4)  =  4  A.  _  (|-J)  Pr 

s  '  VS  f  VS 


5.37 


For  toady  reference,  wo  repent  the  now  simultaneous  sot  below: 

*)  +  s(m)  -  S^(^)  =  -sX  +  ^-dhV 

d?  5.38 

Z-A&hU 


_‘l.  / J£L\  j,  c  /  ^4  \ 

dfc  wj  + 

s(  ^ 


}  4.  /  V9 

'  1  (rtf 


C  77 

■$  M 


®  0 


Immediately  that  If  SmV,  -  0,  which  Implies  %  ' . 
then  the  treatment  reduces  to  a  purely  longitudinal  mode  of  oscillation, 
which  has  already  boon  analyzed  by  Crocco  and  Cheng  In  Ref.  |R. 

lot  us  proceed  with  the  solution  of  Eqs.  5.38.  Rearranging  the 
third  equation  and  substituting  It  Into  the  first: 

ai (>o  / +  ( s  +  it  )( 5w)  “  " 5  ^  +'  ji  " ^  '|- M  +  5  -ji ) 


Atldlng  and  ’Erecting  terms  l„  this  equation  and  also  In  the  second  Equation 
3.  >n  above f  wo  derive: 


-*c  »•- 


25. 


4  /  iJ% 


v  Vo 


7+  S-y)(^+W)  -  s(w*X) 


5,39 


-sV,(v+M-f ) 

i  ( 7*  +  W )  '■  J  ^  +  T  U )  a  "  5  (  V+  '2 ) 

+"  2  S  n  L4 


5.40 


For  con von  I  one o  we  lot: 


D(*) 


(  !$  -  V  +  2abu) 

fe  +  w) 


5.  41 


F  (*| 

E  (*) 

G  (4  ) 


and  then  Eqs*  5.39 


l  S 


V  +  2 

w  -  x 

V  +■  M 


W 


in d  5.40  become: 


+  &fa)  B  *  st  -  sU,  0 


5 1  +  2  S^V\  U 


5,42 


3.43 


Wo  may  eliminate  B>  by  differentiating  Eq»  5.42  and  combining  It  with 
Fq»  5. 33: 

D  a  J}\.  (&l 

where 


5, 44 


Jk(i)  •  si§  + 

di 


Since  the 


s'  F  +  stih 


F-jf  ~i(*+^U] 


light  hand  si  do  of  Eq  •  5.44  is  not  given  explicitly;  we 
will  use  the  method  of  variation  of  parameters  to  solve  this  ordinary  differen¬ 
tial  equation.  Wo  obtain: 

4  2  -CtSwKl&sW,  2  +  5-« 

^  $*’+  S2hh 


£>(*)  =  -  c  i  coaIv  ™ 


■ifi 

'■44? 


29 . 


Wn  can  now  substitute  bad,  for  iron,  Eq.  5.45.  Then  lutegrating 

by  parts  and  combining  terms,  where  the  appropriate  initial  conditions 
obtained  from  Eq,  5.31  i.e.,  f  (o)sO  and  (j-/ o)  A  Z,  'Vr<> 

j{Jl*  — -  t 


or  o 

there  is  obtained: 


ft 


P(fc) 


n 


*  -  Y  +  SU  IJ  =  ■  •C|C0sh\[sl+s^1[n  3  i  5.47 


+■  sf*  E(a')  ttOlMsL+s^h  (i-Z')<h' 

f  £S»  V  |  |  I*—* i.  i.ij  ij..^ 

S’nh  j0  F  (V  )  (-i  “  ?  1 )  c^i?1 

-  s\h  f* G-(i')<5orh{sl+j^  (*-z')<k' 

&h)  [*  O(i0  s;^f?^rh  (w.)  <U' 

Nolimj  that 

l| (*,*')  J3 '  +  -f(«,a) 

wo  may  differentiate  f!q.  5.47  and  substitute  back  into  Eq,  5.42  to  obtain  B 


B(*)  3  “ ^  d*  W  -  pi— 

fs\s*h 


iJ'ft 

.1. 


CiSir)h(?Zrth  &  +  Cj, 


5.40 


(s^sU,  I*  £(*5  ^MsViVh  (a- a')  da 
-  s  (i  F(z')  coihfs^j (a-apofs 


4*  S 


i. 


in  K 


G(V)  Sinb|sl>s*h  (*-*')«!% 


4" 


2  s\>h  j  U(%')  c°at\i|0'+ju, 


We  ran  now  uniquely  determine  the  constants  0|  and  (7 ^  by 
Itoduclng  the  conditions  at  2*  »  0,  From  Eq.  5.31  wo  have: 

Y  fa)  *  W  ( 0  ^  2  Q 

1  V5. 


in 


5.49 


30, 


flnJ  lioncu  « I  mil  I  Innoous  solul  Ion  of  !  qs.  3.47  and  5.40  ,.vl  £  m  Q  yields: 


C,  s 


s4h 


Cu  -  ■fF+s'^ 

and  now  subs  I'  1 1  u  ting  Ihooo  results  Into  r  qs-  13.47  nnd  3.40,  wo  finally 


o!)  I  a  I  n : 


|#‘  «£(v)  sinKfsN-s^w'CwOcI*1 
-3  f*  y  F (*•')  coshiJTVsVW  f*-2r')ol^' 

Jo 


vsx+sV,K  )« 

•t-  2.s%K  (*  i  0  (2*)  cojhfs1>s  Vi  (a-i')nlz 

I  O  • 


*24 

ft 


4^-  P7T7 

»  =  ?  +  sfK  #4;^c^sW 5-  5 

+  »y  -  "4yU  +-  s  [%EfH')c.osb{4Fh(z-z'Jda 
+  '{P+sV  [*  K  F(i’)  sinh(sH7u7  (i-'Adds1 
-  SVi  I*  S  G(i')  Cosli  {s^Fk  (i-Z')tla' 

-2{f;sT;  ( (*  »u(*') 


^  a  sN*  S  nh 


Sjnjo 

5 


Wo  may  also  obtain  a  solution  for  !£&  by  introducing  tq.  3. ‘31  Into  tho 

TO 

third  F,q.  5.30  as  follows: 

-,W  i|  +  |=f;  JVM**>f?5E 

KF(fc')  COiF  fs^t-S^ln 

Q  v 

-  sVk—  1  ff  y<5(«.')  SmbfFfsVW  (i-fc'3cU’ 

-■till  «lll  III!  I  I  III  ■"«■»  I 

i  J  \  J  /O 

-  z  -—2.0*  U  (i1)  coshf^FTTTh  (a-gO'h1 

O  /q 


It  In  c loot  I'hai  Mio  01  dor  of  magnitude  of  tho  individual  terms 
compi  I  sing  l-qs.  5*51,  5.52  rind  6.63  depends  on  tho  ardor  of  magnitude  of  5 
^nd  Ovila  * 

Wo  are  also  intorostod  In  deform  I  n  I  rig  tho  entropy  variation  In 
I  ho  gas,  since  a  knowlodgo  of  this  perturbation  will  bo  necessary  In  apply¬ 
ing  tho  boundary  condition  a  I  the  on  franco  of  tho  exhaust  nozzle.  In  dimen¬ 


sional  form  wo  write: 

i  C  # 


TffdS*  -  db*-  ft  ^p1 

which  may  bo  non-dlmonslonal Izod  to  yield: 

Td5  -  dK  dl 

and  since  dh=  d  T  ,  this  is  equivalent  to: 


!3. 54 


5.55 


T& 


5.56 


V 


and  now  Introducing  tho  equation  of  state  3.8c,  this  becomes; 

S1  =  ?  f  ~  f'-f-  4-57 

In  order  to  separate  the  variables  in  this  equation,  we  take: 


3 


\\ 


e  (*)  4(0  $(9)e 


5.  58 


and  then  Introducing  kgs,  5.1  and  6.58,  (Eq «  5.57  becomes: 

T*  e  i  l  — 

fj  tUM4  »■■"»  «•»  UUMim  V 

“  x  ^  f.  T 

And  now,  rearranging  the  energy  equation,  Eq.  5.29,  we  have: 

iu  Y°  '  mj  *•-« 


?  IV.  % 


4- 

ft 


5. 59 


5.60 
If  - 1 ry5 


This  linear  differential  equation  has  tho  solution; 


e“sWo 


s  4. 

if  yo  C 


dl'  +  coflS-f, 


5.61 


where  the  constant  must  bo  taken  zoro,  since  tho  loft  hand  side  vanishes 


32. 


til'  ^  0.  Rearranging, 

i'  f  *  {-vf  -f  e'1* 


6 .  02 


which  may  bo  I ntograted  by  parts  to  yield: 

.  »J  -L  V  I,  ^"1 


ti  ...  it 


1  _  i  t 
% 


5.6i 


f  =  -tlpVi  ^  +  L-k 

W  %  K  V. 

Combining  17. q s •  5.59  and  5.03  wo  have: 


1C  vt  /, 


e 


d  I  V:v(t')  JJI 


IMXW 


i  (t 


p  ~  =  ~M  f '*  ^  +  ? 

L  K  Vt 


/?  f  ;  ji?l 

„s 

e  ■  •  j 

/o 


Vj  (f) 

Vo 


5.64 


and  now  expand  I  ng  the  form  In  brackets,  wo  obtain  finally, 

"r'f0  “  ~fr  +  =•  ( ft)  ^ 

<•  Vo  \ft)0  -ft  J*. 

I  /£  ,  V  l 

+  vt  I.  e5,  * 


5.66 


6 «  ItllLM^LLDflJB^ 

In  Section  3,  It  was  polntod  out  that  a  rotation  between  I  ho 
burning  rate  and  the  othor  quant  I  lies  would  be  required.  This  relationship 
will  bo  derived  below. 

Wo  recall  that  a  conditioning  process  takes  place  during  the  time 
lag,  and  bonce  Its  duration  Is  a  function  of  the  physico-chemical  processes 
taking  place  during  the  time  Interval  from  tho  Instant  of  Injection  to  I  ho 
Instant  that  combustion  of  tho  given  propellant’  Is  Initiated.  Since  the 
steady  state  combustion  rate  has  an  arbitrary  axial  distribution,  tho  total 
I  lino  lag  will  In  general  bo  different  lor  different  propellant  elements. 

Following  Crocco  and  Cheng,  tho  total  time  lag  will  bo  taken  as  the 


53. 


sum  of  t.i  space  varying  Insensitive  part  Xi  )  and  n  timo  and  space  varying 
sens i I i vc  part  / 

rt(£^0,i)  =  fi(»)  +  Tfa,r  ^t)  "J 


whore  rC  Is  a  function  of  tho  Interaction  index  character  I  nt  Ic  of  the  pro¬ 
pellant  combination.  In  I  ho  nbsonco  of  mixture  ratio  variations,  I  ho  lutes 
may  bo  rotated  to  the  pressure  by  taking 


■f(%/ri  fy’t)  s  f  (^r  (?)  I  +  ^  P 

x  L  p 

where  T  is  the  overall  rate  of  the  conditioning  processes.  Now, 


ducing  Crocco's  definition  of  the  sensitive  time  lag,  wo  have: 


6.2 


Intro 


6. 3 


where  represents  the  quantity  of  energy  requlrod  to  Initiate  burning 

at  station  £*  ,  r  ,  0  at  time  "t  ,  and  the  Intogral  must  bo  evaluated 
following  the  motion  of  tho  particle,  liquation  6.3  may  bo  rewritten  in  the 
following  equivalent  forms. 


where  d  a  denotes  the  particle  path,  and  5 ,  k.  ,  %  defines  the  location 
whore  the  particle  enters  the  sensitive  phase.  We  note  thai  ....cording  to 
the  definition  of  tho  Lagronglen  derivative,  we  may  write,  upon  Introducing 
Fq.  6. I : 


34.. 


t It  (*>'  .Sit) 


\o  (4,nM) 


=  v^o)  +  ^(r)5(0)est 

sjr  - 


rdO  - 
dt 


>N  Wd|(?)ftst 
1  r  ae  £ 


Introducing  Eqs.  6.2  and  6*6  Into  Eqs.  6.4  and  6.5,  wo  have: 


N-91 . 


1  fry 


Vi*  (*')  +  gtf*1)  t  [r'(i')j  <$  [  et4'  ^ 


6.0 


given  by: 

t  (*ir|0) 


/* 


dft" 

Wi'b  +-v;i[i"/t"(z")] 


6.9 


Since  tho  steady  state  solution  is  one-dimensional,  wo  have: 


•P  dt1  »■ 


6.10 


and  furthermore, 
non-un I  form  I l  y, 
becomes: 


since  In  tho  steady  state  there  Is  a  negligible  spatial 
wo  also  have  f  Q  x  ?■:  amsh  Hence  Eq.  6.10 

f  •  t(e)  -  Ecc(i) 


6.  It 


rj. 


and  now  Cq.  6,7  may  bo  written: 


+ 


t  s  \+%  p'[ *'(+'), r'lt'J.eW.t' 


Noting  lha  I-  t(2,r  ,0,t)  =  f  f* )  +  T  '  for,  8,  t  J 
may  nog  loci  second  ordor  I  or  mi.  In  I  q.  6.12  io  obtain: 


T'fc.nM) 


-  ?(S) 


Thin  equation  reflects  the  change  In  the  sensitive  lime  lag  from 
lie  steady  state  value  7(i)  for  a  part Ic I o  which  begins  burning  precisely 
at  time  "i  and  station  it  ,  f  ,  ^  after  having  traveled  through  the 
chamber  with  velocity  [lc‘  (V)  trs(' v‘) ,  $  t1^]  and  having  been 

exposed  to  pressure  perturbations  p1  |  ^'(t1)^  C*  (t*)  f  9  '(t'1),  t 1  “]  . 

Differentiating  Fq.  6.13  with  respect  to  "}:  and  neglecting  higher 
ordor  terms,  there  results: 


/  6,14 

^(sse.t)  =  ~Q-W  : 

dt  p  ( 


or  equivalently, 


6. 1  5 


It  should  bo  emphasized  that  In  writing  those  equations,  wo  am 
dealing  with  a  specific  partlclo  which  enters  Its  sensitive  lime  lag  at 


time  t-Tr  0,i)  (which  coincides  with  station  Ijj  ,  K  >  X )  and  which 
burns  proclsoly  at  time  (the  end  of  tho  sensitive  time  lag)  at  station 


'■>6 , 


2“  )  f  ,9  •  thus,  in  examining  fig,  6,15,  It  Is  cl  oar  that  I  ho 

pressure  perturbation  Is  to  be  evaluated  nt  on  upper  limit  correspond  I ng 

to  I  ho  Irisl  an  I  of  burning,  and  at  a  lower  limll  represent  I  ng  !  he  beginning 

ol  tho  sonsltlvo  lino  lag.  We  will  now  see  why  I  ho  quantity  must  be 

dc 

known  If  we  are  Interested  In  determining  tho  burning  rato  perturbation 

which  can  occur  at  a  glvon  location  In  I  ho  chamber  „ 

•  « 

Let  yvt-  denote  tho  propellant  I  n  Joe  1 1  on  rale,  and  Wb  tho 

propel  I  cm  I  burning  rate,  and  now  consider  the  fraction  of  Injected  pro™ 
po  1 1  ants  burning  In  unsteady  state  between  stations  it  and  r*i  l  ri  it  , 
f  and  rt-cir  ,  9  and  .  Call  this  fraction  and  as¬ 

sume  Hurt  tho  total  time  lag  "[j.  ('4,^ 9, h)  which  this  fraction  experienced  In 
reaching  2,  f  fQ  Is  tho  same  for  all  particles  or  elements  within 
tho  fraction  <fuYi|?(fe||ri 0(t).  Now  this  fraction  was  Injected  at  time 

-  Tfc  as  a  fraction  of  the  Injection  mass  flow  rate 

W\C  .  Then  since  tho  fraction  which  burns  In  time  c^t  was  Injected  In 
I  lino  ift-Tr),  It  Is  clear  that  tho  conservation  of  mass  yields: 


6 , 1 6 


and  hnneo  we  must  have; 


Wi(t-T*  (i,r,0,-t))  { 1-  =  W(qr,0,t) 


0.17 


since  although  7T{  varies  with  location,  It  does  not  vary  with  time.  In 


An* 

tho  steady  state,  — -  vanishes  and  wo  obtain  simply 

dt 


Sfy  l  a  =  covst. 


6.  IQ 


In  an  Investigation  of  Intrinsic  Instability,  tho  Injection  rate  Is  constant 
and  I  lion  wo  also  have: 


£rtit  =  S  YY]i  -  dwb 


6.  Id 


—■•«28RS 


If  wo  introduce  cj:)  or,  Iho  Instantaneous  *  a  I  o  por  unit  volume  al 
which  gas  Is  produced  at  any  point  in  the  chamber,  then 

~  r  cj?  o IQdrdz  6° 20 

by  definition  of  Iho  burning  rate  of  the  fraction  considered.  In  the 
stoady  stale,  this  same  fraction  would  have  burned  botwoon  stations  2 
and  $■  tdi?  ,  F  and  -|-  d F  ,  Q  and  Q  fd&  ,  and  again  by  definition, 
Its  burning  rate  Is: 

(jFdb  (£/?,&)  =  F  5  (% ,r,$)  do  dr  cl%  6#21 


and  now  combining  Eqs.  6.17,  6.19,  6.20  and  6.21  we  find: 

r ^  if)  dr  dll  «  r  $(£,?,§ )dff  J? jl  i*  V  ■ 


In  words,  this  equation  says  that  at  a  given  Instant  of  time  {  , 
the  fraction  of  Injected  propellants  which  burns  at  a  particular  location 
2  ,  r  ,  i f>  Is  a  function  of  that  fraction  of  Injected  propellants  which 

burns  In  steady  state  at  location  2  ,  r  ,  (•?  ,  arid  of  the  variation  In  the 
time  lag  which  occurs  as  a  result  of  the  pressure  fluctuations  during  the 
coordinating  process,  (see  Fig.  6.2) 

Separation  of  the  variables  enables  us  to  write: 

4i*iri i o,t)  =  £60+  &(i)  tYi  $(o)e^  6.23 

whore  ^  s  i('4j) 

and  thus  E q .  6.22  now  becomes: 

r  <^(e)  jfldrdi  +  it(i)  r  t(V)  §(e)  e5*  J9drJi  =  e-24 


r 

da 


r  &L  (iyfil’b)  do  drdU 

oil  di 


3U. 


Lot  us  now  relate  Die  area  el  onion  i  nig  cl  if  which  is  the  steady 
state  area  o lenient,  to  the  area  element  rdOrir  pertinent  to  unsteady  opera 
Mon.  Introducing  the  flrsl  fwo  Eqs.  6.6,  wo  may  write: 


(Y 

|  dr 


r*  V.tr'  [a'.r'fc1),  S'te’U'f*1)  ] 

°  %  [zyM.dWti')] 


where  0  ,  Yq  ,  (9o  urn  the  coordinates  of  the  point  of  Injection.  See 


I g.  6. I  , 


Noting  that  to  £  T  ,  wo  may  neglect  higher  ordor  terms  and 


wr I  to : 


r  ,  F  »  f  *  1  Jt  i  <3.2,; 

;•  M**) 

Replacing  by  the  separation  variables  given  In  Eq.  5. I  we  have: 

r-F  =  /*  Ir(t')  6.27 

,0  %d‘) 

Now  is  actually  given  by  Eq.  6.9,  but  for  the  purposes  of  evaluating 


a  perturbation,  we  may  fakes 


* 


and  hence  Eq.  6.27  becomes:  f  ^  (|^ ^ _ 

r-r  =  ^(r)5(8)e^ffcj)^jL!^Mda'  e-29 

dr'  '• 

on  neglecting  second  order  terms.  On  Integrating  Eg.  6.29  by  parts,  we  ob¬ 
tain,  correct  to  terms  of  0(M), 


r-r 


r  ^  cs  s 

)meiLf'k  '%£>-  ^e-sL  w 

s1  L  r.  % 


Noting  that pQ^pf  and  utilizing  Eqs.  6.6,  wo  rnay  likewise  show  that: 


39, 


r(e-§)  =  H;(r)  dfjjo)^  %  Jk_ 
r  aO  sl 


^r(t)  _  ^  -*(„ 


E'  k' 


ft 


ft 


6.3! 


01  fforont kiting  Eqa.  6.30  and  6.31  wo  have: 

dr  ~  dr 


s*- 


Sl 


6.32 


6. 33 


and  roarranging  Eq.  6.30  givos: 


|p**  y-<* 


*  <k' 


fdr 


6. 34 


Taking  the  product  of  Eqs.  0.32,  6.33  and  6.34  and  neglecting 
higher  order  terms,  wo  find: 

-£%{  f  «■)}(/ 

(  Vq  fto  j 

+  dvfc)  ,  dx 
N/(r-)Jr‘L  i-ij/jrjclr  r 


6.33 


and  then  Introducing  the  right  hand  side  of  Eg.  5.16,  we  finally  obtain: 

I  +  A  SnkY’af  tt£*L  ^  o"S  ^O'kwi/u\otl 

S  l  &  yT  6  W'4'(r)El°)G  |o.36 


J 


Now  Introducing  this  result  back  into  Eq.  0.24  wo  find: 


'•  MlisSS&iteteh!; 


Expanding,  eliminating  higher  order  terms  and  regrouping, 

+ 1 (5)  §  sV'  %  {  %&)  -  tW  i(e)esid 

L  (fo  ) 

-  ( £ )  d  5  ijz  A*  r ,  0,  i ) 

&  ft  ' 


Turning  our  attention,  wo  procood  to  an  examination  of  £?  gjvon  , 

oJi 

hq.  6. 1  ho  From  our  separation  of  1  ho  variables,  wo  may  write; 

?' [**rie»tj  =  Y^mOifeie^  6.:» 

P'[$(vW),k(  )y~(  )M( )]  =?[K  )]t[w  )]$[x(  )]i^r(!l0 


Expanding  each  factor  into  a  Taylor  series,  and  neglecting  products  of 
perturbations,  wo  obtain: 

p1  -  Y(5)t(r)  +K^W  order  , 

•tar^s 


i  t  j 


--  ni®]Wi(9)e>!*'wf 

4*  knher  order  +*ev*nnS 


And  now  substituting  into  Eq.  6.13  wo  have: 
d  Z  u  M  Ot  f  l  D/r\  ( n  F  C  /  r*V1 


si 

6  6.43 


Therefore  fcq.  6.38  becomes: 


41 . 

6.44 


cW  =-  +•  j|(?)(l* 

■h  3k  j  W)  -  g[t(S)]z-iiM}w-W)etiJ* 


) 


Lot  us  now  integrate  each  term  from  zero  to  the  appropriate 
upper  limit,  l.o.  "2  or  H  as  the  case  may  be,  since  there  Is  a  one-to- 
one  correspondence  botween  5*  and  2  over  I  ho  I r  rospoctlvo  paths  of  in¬ 
tegration  when  the  time  Is  held  fixed.  Wo  wrlto: 

lo  ck1  Jo  d?1  Jo  da' 

4  f  (4  $L  ft  (  ?±(v)  -  g-tfiMl  ^(hdz1 

'“s'-  [  ft  &  j 

4  tl'-)i(9)est  |* iL % { i±£)  _  ft[i(aj]e^l  &(*')<!*■' 

ft  1 ,Jt' 


Not  I ng 


that  y%(o)  ™  G  -  ^(o) 


wo  obtain; 


(  W (?• ) 


6,46 


dl1 


From  a  Taylor  expansion; 

V*.(fc)  -  Vi  (^j  + 


+ 


I  •  I  • 


6.47 


and  hence  by  analogy  with  the  ono-cllmonslonal  treatment  In  Rof.  10,  pages 
1 12-1 13,  It  may  bo  shown  that: 


l 


VH§V),?'ln''  +h° IUlt0'  in+°  Kq‘  °'M>  w° flnd  tl"t  +h0  ,actor 


may  bo  o  I  im  Inn  tod.  Ilonco, 


"  ”  *  '”rt‘  -  ™.  «■.,  0,  „ 

'*  -T  *°  - »»>•  is,  ....  „„  *, 

fcrm.,  of  )  or  higher,  tho  remaining  forms  are: 

nj‘  W[,.e-*W-] 

J  dr 

*  7 


-  — „  c/\4  (a) 

s  jr  w 

when  It  Is  assumed  that 


)  ~  ,  —  tLJL,  ,  -L  aZ*. 

y*o  c/a  '  J~ 


'  j£  =  °('j 


In  Section  3,  wo  derived  expressions  for  I f 

y  (ft)  o  6(2;) 

— Tj$-  and  which  aro  given  rocpec lively  by  Eqs.  5.51,  5,52, 

3.53  and  5.65.  Wo  must  thoroforo  examine  the  order  of  magnitude  of  the 
following  Integrals  which  appear  In  the  aforementioned  expressions. 

/o  7, 1 

jj  *  F(a’)  A  (t  -a1)  da' 

X  G(e')  ^  fFTTTw  (*-•£-') dz1 

(W)d*1 

Where  E  ,  F  ,  G“  ,  and  U  are  defined  through  Eqs.  5.32,  3.35  and 
5.41.  Restricting  our  attention  to  the  case  $  ■  0(1),  SnF  *  0  (I), 
an  Involved  term  by  term  analysis  shows  that  for  tho  purposes  of  evaluating 
the  four  Integrals  In  7.1,  iff  ,  tfp  ,  *6  and  VU  may  bo  considerably  re¬ 
duced  to: 


E(fc')  ^  ^  5  -  tf)  Va (a') 


where  ~  %  f 

^  }o  % 

Is  the  first  term  In  Eq.  6.50. 


All  terms  which  do  not  appear  explicitly  above,  yield  contributions  of 
2 

0(M  )  or  higher  after  Integration,  while  those  terms  which  aro  retained 
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yield  contribution,  of  «M).  Furthermore,  5f Y  ,  *U,  *M  and  * 

af 0  a*so  Qac^  0(M)  and  banco  wo  may  wrlfo: 


CCjh  4" ; 


sVs^K  e  + 


f  =  sinKJTvT^rr  ,  , 

v“  S  t 5  "la  2  4-  o(M) 

i  i  | — — — ~. 

<^o  ~  s  CoSM  sN-S^  2:  4“  0(M) 

While  exam  I  no  I' Ion  of  Eq.  5.65  shown  that  tho  entropy  term  |S! 


ft  (&) 
ft 


to  terms  of  0( I ) : 


This  suggests  that  wo  take  as  tho  zeroth  approximation,  correct 


C P(l)  \  (°)  »  r— ■■■■» 

-pT  J  =  cojrnj  sVs^lo  Z 

.‘-ff--  -J^  SiohfssnTo2 

, -  -  j  cosh^j  <rl+r g. 


~s  cosh  { s- 1+ rtK  z 


*Ji 5^  \ ^  -  n 
%  ) 


and  those  may  then  ho  used  ,„  the  evaluation  of  those  terms  and  Integral, 
yielding  contributions  of  0(M). 

Wn  note  that  Lqn.  7.5  constitute  an  exact  solution  for  the  special 
c«o  of  zero  Mach  number,  since  In  that  case,  there  Is  no  combustion  and  the 
equations  coincide  precisely  with  tho  acoustic  solution.  If  wo  rop|aco  the 
exhaust  nozzlo  with  a  closed  end,  so  that  there  Is  „0  outflow,  ihon  *M*J 
rnu.l  he  zero  of  r:  -  ga  =  L  .  The  phenomenon  Is  thus  reduced  to  tho^ 
classical  acoustic  osc, Motion  In  a  cylinder  closed  at  both  ends.  On  sett, no 


a_ 
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tin  d 

o I gonva I uos 


dvolo  /fsN-  s^nh  L  “O 


we  obtain  the 


IjJ  z. 


mx1T' 


1 1 


4-  S  y\\-\ 


(m*  oy  i ,  2. , 3  ■  ■  ■  • )  7.o 


which  are  character I stlc  of  inodes  of  acoustic  oscillations  In  cylindrical 
chambers,  Since  we  have  token  5  ^  ( M)  ,  these  oscillations  must  bo  neutral 
w  i  1  h  wo  1 1  -do f  i  nod  f  i  oquonc  I  os  g  I  von  by  l.q .  7,6. 

If  combustion  tokos  place,  and  the  Mach  number  Is  then  small  but 
finite,  two  modifications  occur,  First,  terms  of  0(M)  must  now  be 
considered  In  evaluating  the  perturbations,  and  further,  the  boundary 
condition  at  £  d  for  neutral  oscillations  Is  no  longer  given  by 

ft 


% 


a  0  ,  but  rather, 


-I-  (ft  7|  -|-  <B  S«h  +■  CL  Ms  s:  O 

7o  7o  To  rc 


7.7 


This  relation  was  derived  by  Crocco  In  Rof,  14,  and  Is  the 
extension  to  throe  dimensional  flows  of  tho  general  solution  obtained  by  him 
I n  Re f ,  5,  T he  coo f f  I c I en t s  Jtf  ©  and  c  arc  complex  fund' Ions  of  the 
frequency,  mode  (  ),  and  tho  nozzle  geometry,  and  tho  perturbations 

at  tho  nozzle  entrance  must  bear  tho  amplitude  and  phase  relationship 
given  by  Eq.  7.7  If  neutral  oscillations  are  to  be  maintained,  Wo  stress 


that,  In  our  case,  wo  could  not  legitimately  take 


To 


h  0  (closed 


end)  as  a  boundary  condition  oven  If  terms  of  0(M)  wore  neglected  In 
Eq,  7.3,  Hence,  since  wo  have  a  new  boundary  Condition  at  z*  *  He, 

It  Is  clear  that  the  values  of  UJ  for  neutral  oscillations  must  now  be 
different  from  the  acoustic  solution,  Eq,  7.6. 

Leaving  LU  for  the  moment  as  the  unknown  eigenvalue  to  be  de¬ 


termined  later,  wo  will  now  set  down  the  expressions  for  the  perturbations 
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which  will  bo  used  in  conjunction  with  bq.  7.7.  I  irst  wo  note  that: 

7.B 
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whore  those  results  aro  obtained  as  a  consequence  of  Eq,  5.6  and  tho  fact 
I luit  combustion  Is  complete  at  2-  ■  'He  so  that 


dV*  (?.e)  »  Afa)  ~  4jiL(%t)  sO 

dtf  ' 


Lot  t 1  rig  S  «  £  U7  (I  n  vest  I  gal  I  on  of  noul  ra  I  1  nstab  1 1 1  ty ) ,  in¬ 
troducing  Eqs.  7.5;  evaluating  the  Integra  1 5  at  an  upper  limit  2  a  E&  ; 

and  making  use  of  Eqs.  7.8;  Eqs.  5.51;  5.5?;  5.53  and  3.65  become: 
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whoro  as  a  consequence?  of  Eq,  7,3 
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the  resell  of  this  first  iteration  is  correct  up  to  tonus  of 
0(M)  provided  llurt  S  a nd  SAi  are  0(1).  In  principle,  one  could  pro- 
cood  with  additional  iterations,  however,  tho  not  rosull  would  bo  to 
iPifroduco  tonne  of  OM*-),  a  ref  I nomoirf  which  Is  not  required  If  tho  Mach 
numbor  Is  sufficiently  small,  In  the  neighborhood  of  0.10,  say. 


Tho  stability  problem  can  ho  stated  ns  follows:  for  a  given 
chamber  geometry,  distribution  of  combustion,  and  exhaust  nozzle,  will 
an  arbitrary  perturbation  of  tho  steady  state  conditions  bo  amplified  or 
damped? 

But,  tho  steady  state  distribution  of  combustion  Is  represented 

by  V*(t)  since  cji  ■  ,  while  the  unsteady  effects  of  combustion 

(Jk 

are  represented  by  the  distribution  of  the  sensitive  and  total  time  lags 
Tfe)  and  Tj.  (%)  ,  and  by  the  interaction  Index  71  ,.  Hence,  stated  con- 
clsoly,  for  a  given  %  ,  T  ,  Tt  ,  (A  ,  ®  ,  C*,  and  71,  ,  will  ar¬ 
bitrary  perturbations  bo  amplified  or  damped? 

Mathematically,  the  answer  fa  this  question  is  given  by  analyzing 
tho  roa I  part  -A.  of  S  #»  .A  i  i  £2  ,  In  practice,  we  need  only 
determine  tho  neutral  condition  under  which  J\~  changes  Its  sign  («A.  »  0  ), 
tho  stability  boundary,,  If  Vj.  ,  Tt  ,  i/fc  ,  and  are  flxod, 
neutral  conditions  will  bo  possible  only  when  a  certain  relation  Involving 
and  tho  Interaction  index  71  is  satisfied,  and  further,  those  neutral 
oscillations  will  'fake  place  with  a  well  determined  frequency. 

Equation  7.7  represents  tho  functional  relationship  between  tho 
three  quantities,  and  slnco  It  Is  complex,  il  corresponds  lo  I  wo  roa I 
equations.  If  for  simplicity  wo  assume  that  the  sensitive  time  lag  Is 
I  ho  same  for  all  elements,  I  hen  wo  have  iwo  relations  between  Ihe  throe 


quantities  i  ,  9b  and  U)  ,  whore  ^  is  Iho  ci  II  leal  valuo  of  I  ho 
sons  It  I  vo  Mine  lag.  This  moans  that  for  a  glvon  value  of  *Y\s  ,  Eq .  7.7 
will  doturmlnu  the  values  of  Iho  time  lag  ^  and  the  frequency  A7  for 
which  neutral  oscillations  can  bo  obtained.  In  other  words,  IEq.  7 ,7  repre¬ 
sents  the  character  I slide  equation  for  the  sol  of  eigenvalue-  £  and 
laj  .  The  most  convenient  procedure,  however,  Is  to  proscribe  the  valuo 
of  U7  and  solve  for  the  o I  gen values  %>  and  &  com pat  I b I o  with 
neutral  oscillations  for  that  value  of  LU  , 

Now  Introducing .Eqs.  7.9,  7.10,  7.11  and  7.12  Into  Eq.  7.7 
and  Introducing  the  notation: 
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where  V)?  and  lo?  are  complex  functions  given  by.; 
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Now,  wo  have  demons  I  ra  led  i  lio  i  lo  -L  and  lo  >  are  uni  quo  iy 
determined  onco  wo  have  specifier]  Hie  rocket  chamber  geometry,  the 
distribution  ol  combustion  arid  the  exhaust  no/zio  goornotry.  Hence 
Eq»  0u2  is  the  final  form  of  the  character  ist Ic  equation  and  may  be 
used  In  the  study  of  the  neutral  stability  of  a  given  rocket  system. 

Lott  I ng, 

Q|  ~  ^3  Re.  -  B,?J 

°2-  "  "“Kim 

wo  may  olimlriato  9X  from  F.q.  0»2  to  obtain; 

f  Ci  z  sinui £  r  Qi  0,6 

Tlio  trivial  solution  buS  «  0,  2TT?  is  discarded,  and 

we  find  that  the  solution  is  given  by  the  simultaneous  sols 

sinw  £  =  =  s>e  A  o.  / 

a^-i-cu1- 

cosual  =  -  coi  A 

+  cuu 

and  hence, 

l  --  d_  [A  +  lK.Tr]  6.8 

Having  determined  one  of  the  eigenvalues,  wo  may  determine  the 
other  by  substituting  into  either  of  the  oqu  I  valent  form's: 

(l  -  CoSut  S  )  - 

^  ~  -  ^  2»  TTrr.  _ _ _ 

^2-Im  (l-ro^cuS)  +•  MnuA 


And  thus  wo  liovo  solved  for  Iho  two  olgonva  I  nos  ?l  and  d 

Wg  note  that'  whon  the  rates  of  the  pliys  Ico-chomlca  I  processor, 
depond  on  the  mixture  ratio,  another  Interaction  Index  exists  arid  a  more 
complicated  form  of  the  characteristic  equation  Is  obtained.  In  the  treat 
ment  of  the  on  I  ropy  wave  instability  analysis,  It  will  bo  soon  that 
(•q.  6.2  Is  a  special  case  of  a  moro  general  relationship  Involving  both 
pressure  correlated  and  mixture  corralatod  offocts. 
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ENTROPY  WAVE  INSTAB  I i  ITY 
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9 .  Tho  Entropy  Wave  Equations 

In  iho  previous  sections,  wo  wore  concerned  with  an  Investiga¬ 
tion  of  transverse  modes  of  combustion  Instability,  in  which  tho  coordina¬ 
ting  mochunlsm  depended  primarily  on  tho  presence  of  transverse  pressure 
waves.  Although  entropy  terms  wore  considered  in  the  I  analysis;  which 
certainly  Influences  the  resulting  magnitude  of  the  interaction  Index  'ft, 
and  tho  critical  value  of  the  sensitive  time  lag  ^  ,  tho  Instability 
Itself  could  not  bo  attributed  to  tho  presence  of  entropy  waves,  because 
these  waves  must  travoi  down  tho  chamber  with  tho  speed  of  the  moan  gas 
motion;  and  hence  a  consideration  of  the  total  period  Involved  indicates 
that  they  can  yield  Instability  only  at  Intermediate  values  of  tho 
chamber  frequency;  i.e,  0(M)  <  uU  <  0(1)  ,  Since  tho  derived  frequencies 
were  of  order  unity  and;  furthermore  could  be  correlated  with  an  acoustic 
mode;  this  would  Imply  that  In  our  previous  analysis;  tho  responsible 
agent  must  be  of  an  acoustic  nature. 

As  wo  pointed  ou+  In  Section  2,  In  which  m  discussed  tho  status 
of  tho  theory;  entropy  waves  may  be  formed  in  either  of  two  ways,  and  If 
both  mechanisms  exist  simultaneous1/;  they  will  reinforce  each  other. 

For  example;  If  we  are  dealing  with  a  monopropellant  motor; 
entropy  waves  can  bo  produced  directly  by  chamber  pressure  oscillations; 
since  neighboring  propellant  elements  wlil  combust  to  a  final  temperature 
determined  primarily  by  the  steady  state  mixture  ratio;  however-;  because 
the  pressure  Is  different;  the  two  sources  of  gas  will  each  have  a  different 
final  entropy.  At  any  Instant  of  time,  iho  disti  Ibutlon  of  excess  pressure 
in  the  chamber  Is  wave  shaped,  and  hence  the  Instantaneous  distribution  of 


entropy  production  In  the  chamber  likewise  forms  u  wave,  and  Hi  Is  wave 
travels  down  I  ho  chamber  and  reflects  pressure  waves  a  I  the  nozzle  exit. 
Thus,  a  closed  loop  now  exists  and  yields  a  mechanism  for  combustion  In¬ 
stability,  since  the  pressure  waves  will  proceed  to  generate  now  entropy 
waves. 

If  wo  are  dealing  with  a  b I  propel  I  an 1  rocket  motor,  mixture 
ratio  variations  may  bo  the  chief  cause  of  enfropy  wave  Instability,  since 
at  those  frequencies,  the  Injection  system  will  respond  to  chamboi  pressure 
oscillations  and  can  produce  off -ratio  mixtures  because  the  oxidizer  and 
propollant  linos  can  respond  differently.  In  this  case,  entropy  waves  will 
bo  produced  because  the  combustion  temperature  will  bo  different  for 
neighboring  sources  of  gas.  Those  entropy  waves  will  travel  down  the 
chamber  with  the  mean  gas  voloclly  and  will  reflect  pressure  waves  at  the 
exhaust  nozzle  which  again  form  a  closed  loop. 

The  total  period  will  depend  on  the  sum  of  the  entropy  wave 
travel  Mmo  downstream  and  the  pressure  wave  travel  time  upstream  plus 
whatever  time  Is  Involved  for  the  propellant  element  to  roach  the  combus¬ 
tion  front.  (See  Fig.  9»|) 

In  the  following  sections,  wo  propose  to  Investigate  combustion 
Instability  In  I  bo  Intermodi  ate  frequency  range  as  It  is  caused  by  entropy 
waves.  Wo  i.all  therefore  consider  the  effect  of  mixture  ratio  variations 
as  might  occur  In  a  liquid  b I  propel  lent  rocket  motor.  For  completeness, 
wo  will  also  consider  the  possibility  of  obtaining  entropy  wave  Instability 
In  the  absence  of  Injector  response.  The  latter  may  bo  termed  Intrinsic 
on  I  ropy  wave  I  nstab  I'll  ty . 

In  Section  3,  wo  derived  a  sol  of  partial  differential  equations 
which  are  applicable  to  a  generalized  study  of  combustion  Instability  in 
liquid  propellant  rockets.  Those  are  repeated  below  for  convenience.  Wo 
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As  wo  have  already  stated,  entropy  wave  Instability  Is  essen¬ 
tia  lly  a  longitudinal  phenomenon  arid  ho  nee  the  following  treatment  will 
bo  ono*’dlinonslonaM  Wo  will  now  lot  our  reference  length  bo  tho  chamber 
length  and  we  take: 


56. 
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5VVV*V>'+  h 

<rvVv«i)  -  -  j  n 

9.2a 

”  T  §f  '  t(bi-^,) 

9r  5a 

P  - 

9  •  4a 

^  di 

9.5a 

“  o 

at 

Hjf  s  =  h/jo  =  h«  +  tl  v/' 

2  * 

9 . 6a 

AM  "0W  W0  prOCOOtl  “*  1)Bfore  "y  Introducing  small  perturbations.  That  Is, 
08Ch  °f  "1a  ‘top0ffldent  varlab!as  “I"  again  be  represented  by  the  sum  ol  a 
Steady  state  space-variable  ami  a  time-dependent  perturbation  so  small 

+l’at  t0mS  hlfjh0r  than  thoss  "lloar  I"  «><>  perturbations  may  be  neglected. 
A  superposed  bar  denotes  steady  state  and  a  prime  denotes  a  perturbation. 
Tlio  steady  stato  and  perturbed  equations  follow  directly. 


d, 

di 

+  (v 


4> 


d  /  _ 


9  *  I  b 


a  ■  = 


It  9-lc 


MaosnJtya 


/  BWI  fcsA 

\  j  —  h  j  j  4*  ~  o 


ir 


Bp 

HIBMH 

Pi 


9,3b 

9,3c 


9  „  4b 


.9,4c 


9 ,  cib 


9.1x1 


9,6b 


9, ,  6c 


Tho  foregoing  sot  of  ordinary  and  per  Mai  d 1 1 f eront la  I  equations 


which  govern  tho  ono-dimonslonal  motion  of  tho  rocket  fluid  system  during 


steady  and  unsteady  operation,  will  bo  util  Mod  In  studying  the  stability 


of  rocket  motors  In  tho  Intermediate  frequency  rango  whore  entropy  waves 


are  operative  In  producing  combustion  Instability. 


10,  Sod* 


the  VarlabJ.es 


In  ordor  to  separate  tho  variables  In  tho  governing  system 


linear  partial  differential  equations,  we  take: 


V'  ‘-1  P(i)e^ 

V  =  (a)est 

f  --  & 

p!  -  $0)est 


p1  = 


(a)  e 

til 


Substitution  Into  kgs.  9.1c  and  9.2c  yields; 


h  s  t 


i(p-vs)- 


,(/) P+\/S)  +■  V’&) 


±  civ 


v/$) 


59. 


Thoso  equations  are  very  similar  to  tho  rosults  that  would  have  boon  ob¬ 
tained  by  sotting  S\oU  «  0  In  Eqs.  *3.26  and  5.27.  Before  substituting 
Into  tho  enorgy  Equation  9.3c;  wo  note  that  It  may  bo  put  Into  the  form 
given  by  l£q .  4.  10  by  Introducing  tho  equation  of  state;  l.e. 


1 1  "t"  ^  i 


M/5VV]  -V  v|p,-Y'-(-()f-i)5VV/ij 

Id  +  hi,  2.  (pv) 


Unlike  tho  treatment  of  tranovorso  waves  In  which  tho  mlxturo  ratio  was 
hold  fixed  and  as  a  consoquonco  was  Identically  zoro,  tho  variation 
of  mlxturo  ratio  Is  now  an  initial  condition  given  by  Eq.  U24,  while  an 
expression  for  Y\^  Is  given  by  (£q.  B29.  Accordingly;  the  enorgy  equa¬ 
tion  may  be  written  In  tho  following  form  upon  tho  Introduction  of  Eqs.10. 
and  B29  Into  Eq.  1 0*4* 


i  [v  jf-TT  +  +  $[-£- Sf  +(x-')fv&j  10,5 
=  £(ft) 

\Ho  note  that  unlike  Eq.  5.29  this  equation  has  a  non-homogonoous  term. 
Wo  will  return  to  the  enorgy  equation  later ;  and  now  lot  us  proceed  to 
the  equation  of  droplet  dynamics.  It  Is: 


V,.  +  (s  +  ^Jk)^  ® 

o  2;  TC 

and  has  tho  Initial  conditions  at  &  *  0  ^ 

(  ^  (o)* 


} 


5 


Equation  10.6  has  Iho  solution 


d£i  r  f  *  fi?.' 

^(t,)  PA^')e("'k)1,  *<l''d4>+™+. 


and  applying  the  boundary  condition 


*](•)  *  *). 


we  may  ol (ml no to 


the  constant  to  obtain:  \  (  ,  , 

w-JL  l0-7 

V; )  ,/u  '  VX(*) 

On  Integrating  by  parts,  the  first  term  yields 

: Ww  »  -  fV^l  +  v£te 

lu  __  s-lvfe  /o  '  di  di 

vfcfr1).  d  $  (f+ik)  f * 1  Jil  7  j3i 

’  fi+Jej1'  S'  j  e  Ji  *»<■■>  j  d'3 

where  we  have  made  use  of  the  condition  £j(o)  =■.  O  •  If  wo  tuko: 


JL  ciJi  <  of  0 
ft  da  -  1 ; 


£  0(M) 

Hl  <  „.... 


then  the  second  contribution  above  becomes  of  O(it')  and  may  bo  neglected 


so  that ; 


Mi)  =  JkMi)  +  n0  W.  £-M)C  fo  io.o 
s 

correct  to  terms  of  0(M) ,  and  where  m  have  not  yet  Introduced  Eq*  13  22* 
The  Initial  conditions  at  the  Injector  face  may  now  be  written* 
*t>( q)  S  /..A  10.10 


*P(q)  *'  ^  “  j:|r.  fo)  *  O  I0J0 

Wo)  *  A  ,  ^(o)"  k/o  iS  Ago  MW’) 

These  conditions  Imply  that  there  is  no  production  of  gas  at  the  injector 
face  and  that  the  liquid  density  cannot  vary  at  the  injector  face. 

The  following  treatment  will  bo  seen  to  parallel  that  given 
for  transverse  waves  except  that  Solo  Is  taken  equal  to  zero.  Com¬ 
bining  continuity  and  energy,  Eqs.  10.2  and  10.5  and  letting 


if  . 


X(t)- 

ifl- 

f)+(iTH)^v| 

_  MM. 

'  WJ 

2 

61.. 

10.11 

y(*)« 

4“ 

~T)| 

\/L^ 

/o 

wo  obtain: 

i  /  i' 
(J-2:  \  % 

>*‘&  - 

-  sX  +  sD^ 

da 

10,12 

On  rearranging  tho  momentum  equation 

10.3  end  taking 

, 

W(a)  =■ 

ii>  V 

i 

%  k 

+  yM 

ft 

10.13 

Z(0  = 

wo  obtain: 

\7  i- 
v  ft 

.  -  < 

£ 

+  y.  'S 
*  % 

d  / 

1  v  \ 

mwiiwrtM*  \ 

,  /i  a 

,  r.  4J 

T  o  “*r“  a 

1-7 

-  S  £1 

- 

10.14 

di  ' 

.ns, ) 

ft 

dOk 

Adding  - 

and  subtracting  terms  In  Eqr>. 

10. 12  and  10.14  wo 

may 

wr  I  to: 

Ui~^)  +  s(*£ +w) 

s  ( W  -X) 

10.16 

fJvVW)  +  s  (fQ 

~  y ) * 

-s(V+Z) 

10.16 

i 

■or  convenience  wo  lot: 

A(V 

( 

~  -  V ) 

%  1  !  ) 

E  - 

- ("  \A/  -  X  J 

10.17 

baf 

vk  +  w)  , 

Ffc) 

=(y+z) 

and  thon  l£qs .  10,1 

1  '3  and  10.16  become 

1  • 

>IA 

•V- 

s'B  =  sE 

10.18 

<JB 

HIP**1* 

t 

s  A  =”sF 

10.10 

dt 

Wo  may  eliminate  by  differentiating  Eq.  10,  IB  and  combining  it 


with  Eq.  10.19: 

(FA 

d?^ 


s’ A  =  Ju 


10.20 


whore 


i  =  s  di  +  slF 


10.21 


Since  the  right  hand  side  of  Eq.  10.20  Is  not  given  explicitly,  use  of 
the  method  of  variation  of  parameters  yields? 

A  (a)  -  -  10.22 

which  may  be  compared  with  Eq.  5.46  after  sotting  «  0  In  the 

latter,  Substituting  for  from  Eq.  10.21,  we  may  Integrate  Eq,  10,22 

by  parts,  Wo  make  use  of  the  condition 


E(i  s  (°)  ~  X f°)  *  JM  + 


■S’  d  S’  Jo 


which  may,  howovor,  bo  absorbed  Into  the  coefficient  of  Jint'iS £ 


and  then  wo  obtains 

A(a)  =  ^  -  V  =  -  C|  cosh -  CiS\\ 

c 

+  s  j0  E(?‘)  co^ksfa-e1)  de 1 
+  s  f F  fi'j  (i -aide 1 


10.23 


Upon  differentiating  Eq.  10,23  and  substituting  back  Into 


Eq .  1 0, 1 8,  there  I s  obto I  nod : 


$  (i)  -  £  +  W  =  C|  c.o.sh sa  +clj\n^sa  10.74 


•s  [4  Ef*1)  c.ojks(-a  -a1)  da 
5  [*  F  (*')  o°fbs  (*  “tdda1 


* 


63. 


Lot  us  determine  l ho  constants  (\  and  Cl  by  Introducing  tho 
Initial  conditions  at  «  0.  From  Eqs.  1 0- 1 1  and  10. 13,  thero  Is 
obtained: 

V(°)  =0  10,213 

\m(o)  =  -  yLv^TM 


and  henc©  simultaneous  solution  of  Eqs.  10.23  and  10.24  at 
yields: 

Cl  ~  o 


X  ■  0 

10.26 


and  now  substituting  those  results  Into  Kqs.  10,23  arid  10.24  w©  finally' 


obtain: 

=  ifY  -  jj  +  to1  7(w)  |  jrinhse  10.27 
+  S  f  *  IT  E(i')  COsUs(i-'i')  <<*' 
t  s  J  *  if  F  M  S'bbi^-fcOda' 


and 


t 

fa 


tfW  +  j  I  +  2~'X fyoty o1  T fu-J)  j  Colitis 


10.20 


ft 


y  E(e')  sibKj(z-e')da 


S  f ^  K  F  (&')  cojhs  (i-z')dz 


Those  two  Integral  equations  will  bo  used  In  tho  solution  of  tho  problem. 
Wo  must  also  obtain  the  equation  for  tho  entropy  variation  In  tho  gas 
before  we  can  write  down  tho  character  1st ic  equation  for  the  chamber. 

From  Eq.  S.b9  wo  havo: 


10,29 


-  1  Wl) 
% 


%) 

fo 


T 


Rowrltlng  the  energy  equation  10.3  we  have: 


Integrating  Eq.  10,30  and  combining  the  result  with  Eq,  10*29, 


so  that  upon  expanding  the  term  In  brackets,  wo  finally  obtain: 


where  for  one-dimensional  flow: 

£  M)  =  6(i)eu 

1 1  *  ItULfilOiM  J3&L£sdM^ 


We  will  now  derive  &  relation  for  the  burning  rate  perturbation 

when  the  mixture  ratio  oscillates.  Following  Crocco  and  Cheng,  the  total 

»*•**■■ 

time  lag  "Pf  's  taken  as  the  sum  of  a  space  varying  Insensitive  part  Xt 
and  a  time  and  space  varying  sensitive  part  X  • 


65. 


Tt(z,i)  =  +  T  ('?.,{•) 


whom  T  Iso  function  of  the  Interaction  Indices  characteristic  of  fho 
propellant  combination.  The  interaction  Indices  may  bo  discussed  In  terms 
of  the  functional  dependence  of  the  factors  controlling  the  rates  of  the 
conditioning  processes.  It  Is  oxpoctod  that  tho  conditioning  processes 
will,  In  tho  cane  of  a  bipropellant  rocket,  depend  to  some  degree  on  tho 
mlxluro  ratio  I**  .  This  may  bo  oxprossod  mathematically  as  follows. 

The  overall  rate  of  tho  processes  at  a  given  location  are  a  function 
‘P(  p;T;  of  pressure,  temperature,  mixture  ratio  and  any  other 

physical  factor  .  We  may  expand  this  function  of  several  variables 
In  a  Taylor  series  about  the  steady  state  operating  condition  where  tho 


local  values  of  the  factors  are  p  ,  T  , 
small  perturbations  P*  ,  T*  #  K  1  and 


and  ^  .  Thus,  applying 
wo  obtain  for  the  now 


procoss  rate: 


f(P(T,r^)  -  ^(fr^T+T'^+K^1)  'n-2 
=  f  (P.T,rp  +  P'  ||  +  T'|f  +  r'g  *  3'  ii 

where  the  barred  quantities  are  to  be  evaluated  at  p  ,  *f  ~T  > 

^  ”•  r  and  vy  e.  pj  .  Lot  us  now  assume  that  tho  temperature  and 


the  physical  factors  are  correlated  to  tho  pressure  and  mixture  ratio, 


T-  T(  p,  K) 

1  an(*  *h0r<>  ^0,*0W3: 
T  -  f  ( i  +  71,  -t-  +  £!  ) 


whom  wo  have  Introduced  I  ho  constants: 


s»f  ,  DT  9-f  ,  QLi  9-f 

®p  9T  3j>  3^ 

El  4.  iX  Sjf  9U  •p.f 

3v*  ar*  st  -Tp  ^ 


66. 


Now  the  definition  of  fhe  sensitive  time  Iqcj  Is: 

*t  *  Tr 

whoro  Ecu  represents  tho  quantity  of  energy  required  to  Initiate  burning 
at  station  £  and  time  "t  .  Tho  Integral  must  be  evaluated  following 
the  motion  of  tho  propellant  element,  however,  since  tho  properties  of 
tho  propellant  element  now  depend  on  tho  mixture  ratio,  we  observe  that 

r** 

td  Is  no  longer  a  constant,  but  varies  with  tho  mixture  ratio  accord¬ 
ing  to: 

c.  B  r.  ( i  4-  /n  „  \^}  \ 

11.6 


Evl  -  Fa 


whoro 


%. 


9  Eg.  f 

r"" 

Ea, 


11.7 


Since  Eg,  Is  associated  with  a  particular  particle,  It  follows 
* 

Ihat  the  perturbation  -is*  must  bo  evaluated  at  the  Instant  of  Injection 


of  that  particle,  Then  Introducing  Eqs.  II » 3  and  11.6  Into  Eq.  11.5, 
and  evaluating  at  time  “t 


rather  It  retains  tho  value  It  had  at  tho  Instant  of  Injection,  It  may 
bo  taken  outside  tho  Integral,  so  that: 

h " 


t  ,n  r 


x-r  t- 


'),t' 


a 


2 IT3E" 


r 


at  M-io 


e  d*“ 


1 


I 


i 
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upon  neglecting  higher  order  terms.  Substituting  back  Into  Iq.  11.9  wo 


find: 


T(e,i)  -  r(%)  =  t(f)  |  ( %  1  ^ jo  [ 

-C  >'  f 


Now  sotting: 


%  = 

flt  =  ^3  “ 


wo  Identify  the  latter  as  tho  Interaction  Indices  and  obtain: 


i?  1  11.13 

|*/o 


This  equation  reflects  the  change  In  the  sensitive  lime  lag 
from  Its  steady  state  value  T^?)  for  a  particle  which  begins  burning 
precisely  at  time  *£  and  station  jj*  ,  when  It  has  boon  Injected  with 
mixture  ratio  perturbation  l*o  at  time  *(;*» )  and  has 

traveled  through  the  chamber  with  velocity  Vfc1  j  anci  l,as  11,00,1 

exposed  to  pressure  perturbat Ions  p'&'hM1]' 

Differentiating  Eq.  11.13  with  respect  to  i  and  neglecting 


higher  order  forms,  there  results: 


a  -  % 


dill* 


Following  the  treatment  In  Section  6,  wo  observe  that  for  one-dlmens lonal 


f low  Eq.  6.1 7  bo comes: 


£yyi|)  ~  <Sv>D  i.  It-Tt  (^)  ] 


(*,£) 

di 


68. 


whore  and  at  o  i ho  fractional  injoctlon  rate  and  fractional 
burning  rate  respectively.  Since  tho  injoctlon  rate  Is  no  longer  a  con¬ 
stant,  Fq.  6.19  does  riot  hold,  and  Instead  m  consider  to  be  tho 


* 

same  goometr  !c  fraction  of 

os 

imi 

Is  of 

Mi  ,  that  Is: 

Srfo  C 

fn  C 

11.16 

Utilizing  Eq. 

6.18,  wo 

have 

$*V\[ 

- 

Sm, 

=  j 

1  + 

1 

]  [  Srhk  ll*l/ 

m  t 

c 

\  &\  L 

1  ) 

and  then  Eq.  11.15  bocomos: 

Srr»t  (*/^)  -  Srvu  (?  )  W/l!iL 


me 


dt 


krt) 


10 


upon  neglecting  higher  order  terms. 

Introducing  ^  as  the  Instantaneous  rate  per  unit  volume  at 

i  ' 

which  gas  Is  produced  In  tho  chamber,  we  may  write 


£mb  (i,i)  -  <t(*rO  dAda  11,19 

by  definition  of  tho  fractional  burning  rate.  In  the  steady  state, 
this  reduces  to: 

Ub  (?)  =  4>(5.)dA4a  ii.2o 


llonce,  combining  Eqs.  11*18,  II. 19  and  11*20  wo  find: 


This  equation  relates  tho  rate  of  burned  gas  generation  at  a  given  in¬ 
stant  of  time  t  and  location  to  the  steady  state  gas  generation 
at  location  B  ,  the  variation  In  tho  time  lag,  and  the  oscillating 
Injection  rate.  Tho  effect  of  the  oscillating  mixture  ratio  Is  of  course 
included  In 

Separation  of  tho  variables  enablos  us  to  write: 

=  $(*■)  +  J|;(Ti)eS^ 


11.22 


69, 


whoro  (&)  ~ 

Thus  Eq.  11.21  hocoitios: 


4 » (£) 


$(?)h  +  -  ji(£)43 


d?.  \  mt  ' 


1 1 ,23 


d5  At 


Before  Integrating  this  equation,  lot  us  oxamlno  tho  expressions  for 

At  U  +  )  /  rtv.  \  ri  .  O"1 .1.  "1  P  4>» t  /«  iii  I  4«U  ITrt  ^  II  I  /!  nn  /I 


il  (&,t )  and  / j5lL)  [i’TlJ  .  Starting  with 

dt  \  VYH  / 

B  24,  wo  obtain  upon  neglecting  higher  ordor  terms; 


Starting  with  Eqs.  ll.Uond 


ji^i)  --  -  &  jr'(S)-  p'CLf-Mjj  "-24 

+  OH^U/W.)  [i-ftfS)]  +  louder*  o^dtfr 

dt  {  p  / o 


Similarly;  wo  may  evaluate  tho  mass  flow  perturbation  as; 

Lt-W)] 


11.25 


Substituting  back  info  Eq.  1 1 ,23  wo  obtain; 


And  now  let  us  Integrate  each  term  from  zero  to  tho  appropriate  upper 
limit;  l.o.  £?  or  £  ,  while  keeping  tho  time  fixed. 


to. 


Wo  may  el Imlnato  the  term  Vf*)-  \/^|)  by  a  Taylor  expansion,  see  Eq.  6.47, 
If  wo  can  first  determine  «  Since  Eq.  11.5  may  also  bo  written! 


11.29 


c 

where  j  Is  the  spatial  location  at  which  an  element  burning  at  •£*•  enters 
Its  sensitive  Hrne  lag,  wo  obtain  for  the  steady  state: 


(a1)  ck' 

%{*') 


11.30 


Into  throe 


Is  a  constant,  and  hence  wo  may  split  the  integral  in  Eq.  11.30 
parts; 


71, 


Upon  neglecting  higher  older  terms,  w©  may  wrlto; 

ft  s  fOt)  -  f  00  (  f*  it'  t  ili  11.32 

;f  \/A(f)  5  '*  Vi.(v)  v£(») 


and  cubst Muting  back  Into  Eq.  11.^1  wo  obtain? 


Now  Eq.  1 1.0  may  bo  rewritten  In  tho  form; 


~  Ea  (  5 )  ^  I  ^  %  3  ^ j  -  %  ( % ,  i)  ■]  "j- 


and  then  combining  Eqs.  11.33  and  11.34  and  rearranging,  wq  find: 


>>  v£(I)  -  '  J^j 

-  j  i  -  |t'j  [t -  Tt M ]  j  |  f  *  ( 1  +  1 1  '3S 


The  term  Involving  may  130  ©  1 1  mlrmtod  by  Introducing  tho 

definition  of  the  Insensitive  time  lags 


On  splitting  the  second  Integral  Into  two  parts,  there  follows; 


1 1 .37 


Subs! i iuting  back  Into  Eq  11.35,  utilizing  the  do  f  I  it  it  Ions  Cq.  11.12,  and 
combining  forms,  wo  finally  tier i v«: 


72. 


This  equation  yields  the  shift  In  the  bur  nine]  station  from  i  ho  steady  stato 

K«3  I 

location  if  for  a  particular  propellant  trlomoni  Injected  at  time  t^X-L 
There  are  three  contributions  to  this  shift.  The  first  term  on  the  right 
hand  side  gives  the  effect  of  a  change  In  droplet  velocity,  while  the  second 
and  third  terms  correspond  to  a  change  In  time  lag  due  to  pressure  oscl  I  ja-». 
t Ions  and  mixture  ratio  perturbation. 


introducing  Eq,  11.38  info  Eg.  6J.7  there  results; 


arid  now  substituting  this  result  Into  Eq»  11.28,  we  obtain  the  form; 


On  examining  this  result,  wo  observe  that  the  perturbation  In  the  gas  flow 
consists  of  tho  contributions  of  two  groups,  each  consisting  of  throo  forms. 


73. 

T  ho  first  .jf'ou  p  of  terms  may  be  cal  loci  tho  t 1  mo -w  Iso  coni  r  I  but  I  on,  and  re¬ 
presents  trio  of  foot  of  tho  perturbations  In  tho  local  burning  rato,  which  In 
turn  Is  duo  to  tho  variation  In  tho  'Mine?  lag  duo  to  tho  pressure  oscilla¬ 
tions,  tho  Initial  perturbation  In  mass  flow  and  tho  time?  rato  of  variation 
of  the  mixture  ratio.  Tho  second  group  of  forms,  which  may  bo  callod  tho 
space-wise  contribution,  represents  the  effect  of  the  displacement  of  tho 
location  whore  a  given  propol lant  element  burns. 


And  now  using  Eqs.  10. I,  10.9  and  the  results  of  Appendix  0,  and 
then  dividing  by  (fi  ,  wo  have: 


Those  terms  may  he  simplified  by  considering  an  order  of  magnitude  analysis 
of  the  six  parts  which  contribute  lo  the  source  term.  Tho  maximum  local 
value  of  Is  assumed  to  bo  0(1),  and  In  addition,  It  Is  assumed  that 


£ 


iJ 


Jv* 


ft  ■  f.  '  %  '  % 


%  d  « 


*  0(0 
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Wo 


also  nolo  that  If  Gr(uj)  >  H  (wi)  ,  J(uj)  and  /V|(ia))  are  each  0(1), 

W  I  |/\A  l  /  I  .1 


then  Appendix  R  shows  that  ~~r 

r* 

order  of  a  perturbation,  which 


yyv  \l  •  i 

— and  tA/i  are  each  of  the 

mi  %n  ^ 

In  agreement  with  our  understanding  of 


these  terms. 
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S  I  nee 

matt or  how  large 
wo  may  show  I  hot  t  ho  first'  lonn  In  Eq. 


must  Integrate  lo  o  quantity  of  0(M),  no 


r  g(*i«h 

^(l‘)  I',  locally,  on  taking  a  Taylor  expansion  In  V I -l)  , 

ck  V?, 

41  bocomos; 


1 1 .43 


whore  the  Integral  yiolds  a  term  of  Q(M) .  The  second  torm  In  Eq,  11.41  Is: 


and  Integrates  to  0 ( M)  and  cannot  bo  simplified,  Noting  that: 


the  third  contribution  to  tho  sourco  torm  bocomos: 


which  Is  of  Q(M) .  Tho  fourth  terms  Is: 

V,(0  £VvLr(i)GM  =  °(^) 

*  (.la 


On  Integrating  tho  fifth  torm  by  parts  we  obtain. 


whore  the  bar  represents  a  proper  mean  value  so  that  when  S  s.  0(1),  the 
fifth  term  becomes: 


s*-  Vo  At 


And  now  taklnn  I  he  last  form,,  wo  have? 


7*3 . 


.Hid  dll  terms  of  O(M^)  or  h I ghor ,  have  boon  nogloctod.  Wo  conclude  that' 


under  the  present  assumption! 
of  0(1). 


3,  W-  Is  of  0(  I)  local  ly,  if 

r<  l" 


1 2 -  Solution  by  Iteration 


In  Section  10,  we  derived  expressions  for  (£^r)  ,  and 

>  which  are  given  respectively  by  Eqs.  10.27,  10.28  and  10.32. 


I  9 

We  must  thoroforo  examine  the  order  of  magnitude  of  the  following  Integrals 
which  appear  in  the  aforementioned  equations: 


12.1 


where  £(i')  and  |r(v!,)are  defined  through  Eqs,  10,11,  10.13  and  10,17, 
Restricting  our  attention  to  the  case  $  «  0(1),  an  analysis  of  the  contri¬ 
buting  terms  shows  that  for  the  purposes  of  evaluating  those  Integrals  to 
within  forms  of  0(M),  we  may  take: 


I  U" 


.3 


\  ( i.  w't 

*£(»)  =  0-*)v(*')  +  L^e-S>»  p^j  Jv^.) 

\  r  —  (ji.1  12.2 

x  F(i')  =  x%  QM  +-  Jrjj  |  h(w) -nlM  sG(w)j  £ 

-M  JMe'5(ip.>  j4f 

whore  j  <-*£ 


(*  J5mw 
je  ft 


and  Q(fc')  has  boon  defined  In  Eq.  11.4  7,  Wo  nolo  I  hot  both  FM  ond  Y (v) 

(5  /al\ 

which  depend  on  ore  locally  of  0(1)  because  of  t ho  contribution  of  the 

term  given  by  Eq,  11.4b,  but  upon  Inlograliun  this  term  goes  to  0(M)  since 
Is  0(M),  and  consequently  ffc1)  end  Y(Z‘)  must  demonstrate  the 
same  behavior,  Thus  all  horms  which  do  nol‘  appear  explicitly  In  Eqs,  12,2 
and  12,3  yield  contributions  of  O(M^)  or  higher  after  Integration,  while 
those  which  are  retained  yield  contributions  of  0(M). 

Since  VgQ  a  0(M),  the  product  ,fi  0^),  dnfJ 

Eqs,  10,27  and  10.28  reduce  to; 

=  xYt*)  -swkss  +  s  j[  [if  EH  tojfu (i-i') 

+  s  J*  [ if  F  H  ( a 

v W(»)  +  cosKss  -sj^  f(i,)  "mV> W s  (*-*•) 1*1 


% 


)  mks  (fs-z')'fda1 


12,5 


2.6 


where  from  Eqs.  10.11  and  10.13,  wo  have  correct  to  0 ( M ) ; 

*V( H)=  -  (fV(e)^y  +o(Mx)  »  o(') 

Jo  Vo 

yW(«)  =  +  Q|(^i)  a  o(m) 

Y>o 


12.  7 


n. 


And  now  s  i  nee  Siv^Vv.>rf.  «  i  SwauJ  £ 
row i  i to  I  q s .  12.3  and  12.6  in  ,'ho  form; 


o-' 


d  (  oilo  Stl:  “  cojuOj?  wo  may 


yj/(z  )  -  4  0(M)  =  tJimvZ  4  o(m)  12.8 

^  V2o  vy*0 


5-  (t )  -  coslosz-  4  o(m) 

¥v> 


Cojuo-2  4  O(M) 


9  0 

(n.  f  2 


will  I o  examination  of  Fq .  10.32  shows  that  the  on t ropy  form 

*«£) 

% 


o(0 
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This  form  of  t ho  aquations  suggests  the  following  Iteration  procedure  . 
First  neglect  the  forms  of  0(M)  with  respect  to  forms  of  0(1),  and  then 
utilize  the  zr-*oth  order  solution  to  evaluate  the  higher  order  terms.  The 
resulting  zeroth  ordor  solution  Is* 


/ir£(i)\(o) 

%  J 

V(0U° J  . 


i  sinoui* 


12, 


C  o  J  (At  % 


correct  to  terms  of  0(1),  where  from  Eq„  11.46  the  form  which  Is  local  I /  of 
0(1)  Is: 

'*  ck'  (i  «  p'  ,-k" 

0-JrtJ-o  fcfen) 

12.12 


(°) 


%M  f 


yv) 


'o  v/(*0 


As  we  pointed  cut  In  Section  7,  when  the  Mach  number  Is  Identically  zero, 
there  Is  no  combustion,  -..El? 

»  0  ,  and  there  Is  no  outflow,  so  that 

% 

Fqs,  12.11  coincide  precisely  with  the  acoustic  solution.  That  Is,  the 


phenomenon  Is  reduced  to  the  ono-d I  mens  Iona  I  oscillation  In  a  cylinder 
closed  at  both  ends.  At  the  chamber  exit,  9;  =>  I,,  so  that  the  corresponding 
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I  gen  v©  l  ues  sat  I  s$y  flte  equal  ion ; 


UJ  s  MTT 


12. 13 


character  I  s+  Ic  of  organ-pipe  oscillations.  t  nrlhonnoro,  the  one  I  llal  Ions 
mus'l  bo  neutral  at  those  we  1 1  -def  I  nod  frequencies. 

If  combustion  occurs,  and  'the  Mach  number  is  small,  but  different 
from  zero,  terms  of  order  M  ore  added  Into  Eqs.  12.8  arid  12.9,  but  of 
greater  significance  Is  the  fact  that  the  boundary  condition  at  i?:  ■  t, 

in  t\r\  I  AKtrirtK"  r-i  I  .  j..n  K . .  ^  ^  ^ '  A  l\  i  1 4  n I  U  .  .  n  I  .  . . 


s  no  longer  given  by 


«•  0,  but  rather  by; 


a.  04\Ve  It  +  (lA  Ve 
ft  ^  %  I  % 


12.14 


This  relation  Is  the  one -dimensional  form  of  Eq»  7.7  derived  by 
Crocco  in  Ref.  19,  (see  Appendix  A).  Wo  note  that  f>(^  and  /3va  are  com¬ 


plex  functions  of  the  frequency  and  nozzle  geometry,  and  hence  It  follows 
that  since  w©  have  changed  the  boundary  condition,  the  values  of  caJ  for 
noutra I  oscillations  are  no  longer  given  by  £q«  12.13,  even  If  terms  of 
0(H)  are  neglected.  Leaving  uj  for  the  moment  as  an  unknown  eigenvalue 
to  fee  determined  later,  m  may  now  set  down  the  equations  for  the  perturbs 
tlons  which  will  bo  utilized  In  conjunction  with  |q.  12.1 4.  First  wo  note 


that; 


V  V  (&c) 


% 


if  \  (*) 

vsj 
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+■  t [  d v (ihte 


h  V(0  C .0,s  uj 


i  W  (*«.)  =  x\  i  XPt  \  CO  E  _  Zl  \jt  VJ 

'  ft  / 
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whore  those  results  are  obtained  from  I  gs.  12.7  and  11,46,  and  the  fact 

that'  combustion  is  complete  at  £  «  a  1/  so  that  ]  0. 

ciZ  /  Ze. 

Evaluating  I  ho  integrals  at  the  chamber  ox  It,  wo  obtain  from 
Ec|S»  12,3,  12,6  and  10.32: 

=  k y (t V^e)  - ti'muj  4  [ir V*1) cor u^i - 1  12,16 

\  ¥0  f 

-iaJ  I  [)f  F(°Va')  Simv(l-S')]^' 


(  ^  ^  )  +  CoSl*J  +  W  I  [*E(V*')jlMW(l-H')]c(i'  12.17 

-  tw  cosi^fi-i1)]^1 


p7  yv  =  12. 


+ 


■Wi 

Ve 


>0 


11 '  '•  V(ii  y  (*')  w  SiVu/J«'c#8 


10 


4* 


Ve 


e"lwio^;)elW  |  da. 


4y  (i')cid 


Is  noted  that  r  J  and  f  “7 2T"  )  are  given  by  Eq*  12.11 
and  are  to  bo  utilized  In  evaluating  tho  Integrals  of  icE  frj  and  *F(>) 


and 


wherever  appl  i cab  1 0 . 


3.  D I scussi 


We  may  now  stato  the  stability  problem  as  follows:  for  a  given 
Injector,  chamber  geometry,  distribution  of  combustion  and  exhaust  nozzle, 
will  an  arbitrary  perturbation  of  I  ho  steady  state  conditions  bo  amplified 
or  damped? 


Note  that  comparison  with  Section  0  shows  that  several  now  factors 
ha vo  boon  Introduced  Into  tho  mathematical  formulation#  That  Is,  I  ho  steady 
state  dlsti  I  but  Ion  of  combustion  is  still  determ  I  nod  by  (i.)  ,  but  now  in 
considering  tho  unsteady  effects,  duo  regard  must  be  taken  of  tho  Injection 
system  geometry  and  character  I  st  Ic  time  (X-  ,  and  also  tho  effect  of  the  second 
Interaction  Index  7fl  .  The  revised  statement  of  tho  stability  problem  there¬ 
fore  becomes;  for  given  &  ,  Vf'i)  ,  T(^j  ,  X{(l),  ^  and  71  , 

will  arbitrary  perturbations  bo  amplified  or  damped?  As  before,  wo  will 
determine  the  stability  boundary,  Lo.,  wo  will  solve  for  tho  neutral  condi¬ 
tion  J[,  0,  5  «  Of  . 

It  0<  ,Vfe)  )  Tt  ft)  t  and  ^>v\  are  hold  fixed,  those  neutral 
conditions  aro  possible  only  when  a  certain  relat  ion  Involving  Tf*0  and  the 
Indices  and  %  Is  satisfied,  and  they  will  take  place  with  o  woM  defer- 
mined  frequency  UJ  .  Eq«  12,14  represents  tho  functional  relationship  between 
tho  four  quantities,  and  slnco  It  Is  a  complex  equation,  It  represents  two 
real  equations  relating  tho  sensitive  time  lag  Z'(l)  and  neutral  frequency  (A) 
to  tho  Indices  Tit  and  %  .  For  simplicity,  we  will  suppose  that  tho  sensitive 
time  lag  Is  tho  same  for  all  propellent  elements,  whore  <5  represents  the 
critical  value  of  the  sensitive  time  lag.  Thus  for  given  values  of  7/i  and  71  , 

• wm 

It  Is  possible  to  determine  tho  values  of  the  time  lag  £  and  the  frequency  UJ 
for  which  neutral  conditions  can  bo  obtained,  Eq.  12.14  Is  the  characteristic 
equation  for  tho  set  of  eigenvalues  and  UJ  .  However,  as  Indicated  In 
Section  8,  the  most  convenient  procedure  Is  to  proscribe  the  value  of  UJ  , 

M-" 

and  then  to  look  for  the  eigenvalues  %■> ,  %  ,  and  £  compatible  with 
neutral  oscillations  far  that  value  of  uJ  . 

Now  substituting  Eqs.  12.16,  12.17  and  12.18  Into  Eq.  12.14,  we 
ob to  In,  ( see  a ppond lx  C) s 

»j 

h,  %l  *  h^(i-e-Lwf)  +  ' 1 3 


13.1 
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where  V\\  ,  1^7,  mid  lo  3  are  comp  I  ox  numbers  which  aro  a  un  I  quo  function 
of  f ho  frequency  onco  the  rockot  chamber  geometry,  the  Injection  and  exhaust 
systems  and  the  distribution  of  combustion  have  boon  specified.  Therefore 
Eq.  13.1  Is  the  final  form  of  the  characteristic  equation  of  *1  ho  rocket 
chamber  and  may  bo  used  to  Investigate  the  stability  of  a  given  rocket  system. 
Note  I  hat  this  equation  Is  somewhat  more  complicated  than  Eq.  8.2. 

Separating  tho  real  and  Imaginary  parts,  Eq.  13.1  becomes: 


QVl  S  h  |  Re  : 
wS  Kim 


%[h,.  Cojw?)  +  $w\uj1  "]  4  l^3i^ 


Essentially  then,  wo  have  two  simultaneous  equations  of  the  form: 


Tj  ( $,  m)  so 

h  w)  "  0 


13.3 


Lot  us  first  eliminate  9h>  ,  and  then  we  obtain  tho  single  equation: 


% 


13.4 


mm  — — 

ku(i~CoSu;i‘ )  -  bjSinwJ 


where  bi ,  K  and  are  roal  numbers  and  are  constant  for  a  given  value 
of  tho  frequency, 

b|  =  b|R<.h3,im  -  biXwfi3Rc  15,15 

b  j.  b  I  Xw\  b  x  "**  U|Re  1^3  Z.  X100 

55  Iw*  Via  I»n  “4  1-8  l  Re  Va % R.& 

Wo  can  Investigate  the  behavior  of  9t  In  tho  Uj£  plane.  It  Is  Immediately 

KM 

nolod  that  %  must  bo  a  periodic  function  of  l/j£  .  By  sotting  tho  denomlna- 

*a» 

tor  of  Eq.  13.4  oquol  to  zero,  we  may  determine  tho  values  of  wi  for  which 
%  becomes  Infinite.  Wo  have: 

\>’1  (|-  CoXuJ  s)  ^  $W8bJ  <5 


13.6 
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Since  has  a  por  ladlclty  of  '  ll  In  I  ho  IAJ  S  piano,  wo  may.  confine  our 


as  wo  I  I  as: 


whoro 


Iho  region 

o  §  w  S  $  2tr  . 

Solution  of  [1 

,q6  1 3.6  y lei d 

UJ  S  ~ 

2  K¥ 

((<=0,  1 

•  13.6 

OJ  <£  ~ 

■  B 

1 3.  7 

cos  B  = 

b3  - 

bj^  t 

13.8 

sin  B  = 

2  k*  bj 

W-+  ^ 

Since  t  —  approaches  zero  faster  than  Sinu^<(  ,  7L  will  tend  to 

Infinity  with  the  sign  of  -b'/b,  as  U/S  approaches  zero  from  the  right,  and 
%  will  tend  to  Infinity  with  opposite  sign  ns  ui£  approaches  2TY  from 
the  left.  Hence  U)S  ■  0,  and  uij  ■  2.T  constitute  vertical  asymptotes  for 
fl  .  The  third  vertical  asymptote  of  yt  Is  given  by  Eqs.  1 3. 7  and  13.8. 
Examination  shows  that  must  approach  ini  I n II v  with  the  same  sign  as  at 
the  origin  to  the  left  of  this  asymptote  and  with  opposite  sign  to  the  right. 

i 

If  ot,  Is  a  continuous  function,  It  must  take  on  a  stationary  value  between 
each  pair  of  asymptotic  values.  The  two  values  of  i/jS  for  which  9l>  Is  sta- 
tlonary  may  be  determined  by  setting  r-g.  *  0.  Accordingly, 

d(t4 ) 

b|  Pi  -  r“l  13.9 

-- — — ■ — — - - - -pjCoSwiJ  -  O 

C^i(l'"Co^^)  “  ba  SiowjJ  J  1 

and  honco  assuming  that  Is  not  Identically  zero,  and  that  the  decomlna- 


tor  does  not  vanish,  wo  find  that  '/L  Is  stationary  when 

***** 

bliinwi  =.  b^CojiAJcf 
There  are  two  solutions  to  this  equation. 

f  ~ - — - - 

}  _  + 

(sinwS|  =  ka _ 

\  WL  +  h2- 
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1 3. 1  I 


03, 


COS  <aJ  o  ^ 


b  i.  -I"  b  / 


S  ivy  u)  S 


it  follows  that  4,  and  uii^  will  No  In  the  first  and  third  quadrants, 
oi  In  I  ho  second  and  fourth,  rospoct  I  ve  I  y .  Thoso  results  are  pictured  in 
Figure  13.1  and  13.2. 

It  Is  c  I  oar  that  if  we  wish  to  obtain  tho  minimum  value  of  'fL 
compatible  with  neutral  oscillations  at  a  given  frequency  ,  we  may  take 

either  or  as  given  by  £q.  13.11  or  13.12  and  substitute  back  Into 
kq.  13.4.  Only  r ho  positive  values  of dl have  physical  significance. 

With  these  values  of  %  and  S  ,  we  can  also  substitute  Into 
Eq.  13.2  and  obtain^.  Thus  OH,  Olm in  and  &  may  bo  determined  as  a 
function  of  W. 

Lot  us  examine  tho  function  'ftiftAjfi)  more  closely.  We  may  pro¬ 
ceed  by  eliminating  01  f  rom  If q s .  13,2  tin d  wo  obtain: 

jj  ^  1?4 ( t  “  COSujS  )  -  bs  S\0u4 


where 


bL(l-C0S vd<?)  -  b^ 

H  zRe  h  3  Re  +  hi  Im  h  3  3>n 


13.13 


13.  M 


since  the  denominator  of  Eq.  13.13  Is  either  finite  or  zero, 
can  vanish  only  when  the  numerator  vanishes.  Tho  numerator  vanishes 
when  4  w  0  or  2Tf  ,  but  the  denominator  likewise  vahlshos  there. 
However  wo  have  seen  that  will  vanish  faster  than  $,y>i4  ,  and 


hence 


uJ  &  "*■  XTT 


13. 115 


84, 


But  since  c)  «  0  ,  at  »  0  ,  wn  find  "  C>Q  w  I  I  h  §  I  gn  ^  ,  while 

<*  wS*  w.(it)  =  lif tj  •  _ 

The  numerator  also  vanishes  whon  UJ<S  la  given  by: 

(’  2  ~  b^?' 

C05i/J4  “  - ; 

b/|  *  T~  *“  13.16 

SinujJ  - 

3“  -r  Ks2, 

Since  thn  denominator  of  Eq*  13,13  doesn't  vanish  at  this  value  of  yJ S  ,  In 


general,  and  since  o  Is  finite,^  vanishes  at  most  only  once  at  the  value  o< 
U)S  given  by  Eg.  13.16. 


.  i  r  :*vSn*ff  \J 


85 . 


IV,  results  and  conclusions 


Slnco  wo  aro  dealing  with  a  linearized  analysis,  all  modus  of 
oscillation,  Including  the  standing  wave  and  traveling  wave  forms,  can 
exist  simultaneously  and  Independently.  However,  all  modes  have  a  unique 
frequency  for  a  fixed  geometry  and  hence  may  be  Investigated  Independently 
of  each  other,  Tho  computation  of  tho  stability  limits  for  a  particular 
rodkot  motor  proceeds  directly  once  tho  chamber  and  exhaust  nozzle  geometry, 
Injection  system  and  the  steady  state  distribution  of  combustion  have  boon 
described,  as  discussed  in  Sections  Q  and  13. 

.  Thus,  by  way  of  illustration,  a  typical  procedure  would  be  as 

follows:  specify 

(a)  tho  chamber  geometry,  length  and  diameter, 

(b)  tho  Injection  system, 

(c)  the  subsonic  portion  of  tho  exhaust  nozzle,  arid 

(d)  the  steady  state  distribution  of  gas  velocity  In  the  chamber. 

In  connection  wlJ'  The  above  items,  It  Is  noted  that  the  drag  coefficient  of 
the  droplets  ,  will  depend  on  the  droplet  diameters  and  hence  on  the 
Injection  system,  Item  b,  as  well  so  fho  viscosity  of  tho  gas.  It  Is 
sufficient  to  specify  only  tho  subsonic  portion  of  tho  DoLaval  nozzle,  Item  c, 
since  the  flared  supersonic  portion  has  no  effect  on  the  chamber  oscillations. 
And  finally,  wo  note  that  since  tho  description  of  tho  combustion  process  Is 
obtained  equally  well  by  proscribing  either  the  burning  rate,  or  tho  velocity 
distribution,  for  convenience  wo  may  proscribe  iho  latter  taking  caro  lhat 
the  axial  component  of  tho  steady  state  gas  velocity  Is  zero  at  tho  Injector 
end,  arid  has  a  vanishing  spatial  derivative  at  the  exit  end  of  the  chamber. 
This  last  condition  must  bo  met  If  ono  wishes  to  satisfy  Jho  requirement  that 
combustion  Is  complete  within  the  chamber. 


Bo, 


Returning  to  I  torn  b  above,  It  Is  observed  thal  many  different  types 
of  Injector  response  are  Included  in  the  analysis  of  Appendix  B,  Including 
for  example: 

(I)  cav I  tat  lug  venturi  Injectors, 

(II)  matched  Impedance  Injectors,  and 

(III)  mismatched  Impedance  Injectors. 

An  Ideal  eavltatlng  venturi  Injection  system  would,  of  course,  have  zero 
response  to  chamber*  pressure  oscillations,  and  hence  result  In  the  following 
s Imp  I  I  f  leaf  Ions: 

C4(uj)C)Vi  =  HMcV.  =  J  (^)c.v,  =  M(i*j)c,v,  *  O  14,1 


If  wo  define  a  matched  Impedance  Injector  as  one  for  which  tho  mixture 
ratio  alone  does  not  oscillate  (minimized  entropy  wave  effects)  when  tho 
chamber  pressure  oscillates,  then  the  analysis  of  Appendix  B  shows  that  the 
two  conditions  which  mu  t  be  fulfilled  are: 


14,2 


If  only  tho  first  condition  In  14,2  Is  satisfied,  then  the  oxidizer  and 
fuel  lines  will  havt  the  tamo  phase  lag,  but  the  Injector  wlli  nevertheless 
produce  an  oalcl latlng  mixture  ratio.  This  result  follows  because  the  twin 
conditions  of  equal  phase  lags  and  equal  amplitudes  must  bo  met. 


Tho 

relations 

h  1  pj 

?  in  Eqs 

,  14,?  may 

also  be  expressed 

P  Rox 

=  Pkp 

f° 

.  ^ 

' 

f° 

(i  tt. 

/\(a)  - 

/  A  \ 

1  f 

1 

"  y  r 

A(»)  j 

I 

oX 

(n  oh ox 

K* 

(Ao 

1 

(Ao  ) 

14.3 
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wharo  as  before,  Iho  subscripts  ox  and  f  refer  to  I  ho  oxidizer  and  fuel 
respectively.  If  all  three  conditions  are  satisfied,  then; 

Gr(^J)/v\,|.  ~  M  (iAj)  (v\,i .  -  0  14,4 

to  within  terms  of  0(M)  . 

Thus,  for  cases  i  and  II,  Eq.  C  21  shows  that  hj  Is  Identically  zero 
and  hence  this  moans  that  when  the  mixture  '-atlo  Is  constant,  the  character  1st Ic 
Equation  13.  I  can  ho  reduced  to  the  same  form  as  Eq.  8.2  ovon  when  there  Is  some 
other  form  of  coupling  between  the  chambor  and  the  food  1 1  nos.  Case  III,  of 
course,  corresponds  to  the  general  case  of  a  b  I  propel  larit  Injection  system 
with  arbitrary  relative  phase  lag  between  the  oxidizer  and  fuel  line  and 
arbitrary  amplitude  response  for  each  lino. 

II  Is  noted  that  for  sufficiently  high  chamber  frequencies,  the 
functions  (Jfuj),  H  (w),  J(wj)  and  Ml^)  are  negligible  and  h|  again  vanishes. 

Thus,  to  summarize,  the  characteristic  equation  of  the  chamber  reduces  to 
the  form  given  by  Eq.  6,2  when  the  mixture  ratio  Is  constant.  This  condition 
exists  generally  In  a  monopropol lent  motor,  and  In  a  b I  propel  lent  motor  with 
an  Injector  corresponding  to  cases  i  or  II,  and  also  as  Just  observed,  In  a 
chamber  In  which  only  high  frequency  oscillations  are  present. 

Once  the  propellant  Injection  velocity  Is  utilized  In  solving 
Eq.  4.11  for  ,  Eq.  3.7  may  be  rearranged  to  give: 

%)  -  Ve  •"  VO*)  |4(5 

since  \jv,B  \/ ,  .  It  Is  then  possible  to  evaluate  all  Iho  Integrals 

leading  to  the  determination  of  the  three  complex  quantities  I'M  ,  n ^  and  lq 
as  a  function  of  chamber  frequency,  arid  then  to  determine  the  stability  limit* 
of  the  chamber  as  given  by  t ho  eigenvalues  and  X  .  In  closure,  It  Is 

remarked  that  the  complexity  of  the  Integrals  used  In  evaluating  I'd,  Iq^  and  K3 


on. 


procludos  any  but  numerical  evaluation  of  those  functions. 

Lei  us  now  consider  t ho  two  rocket  motors  dof inod  below: 


Chamber  No. 

Typo 

re* 

L* 

Zc 

4 

'  4h 

<! 

r-* 

10 

JL%\ t 

1 

Short 

2 . 000" 

2.000" 

1 . 000 

0.846" 

0.  1 

0.720" 

2 

Slier  t 

2,000" 

2.000" 

1.000 

1 .  1 50" 

0.2 

0.910" 

whore  Is  obtained  from  Pq,  A  9  for  !<  »  1.00  and  where  I  he  steady  state 

distribution  of  chamber  velocity  Is  given  by  : 

o  <;  £  <:  Q.Z  V H  ~  O 


o,2t(  s  2  4  o.s  z.e  %-  — fa.-o.2ae) 

o.3  Z-e 

•Wit 

O.'jZc  £  Z  *  He  VH  - 

and  whero/$  13  0,15  for  Vp  «  0.013  ( See  Figs*  4.2  and  14.1). 

Since  wo  have  already  out  1 1  nod  all  of  the  fundamental  concepts,  solved 
for  the  eigenvalues  In  Impllcll  form  and  briefly  discussed  the  procedure  to  be 
followed  in  evaluating  all  of  the  component  Integrals,  wo  need  merely  remark  Hint 
we  will  Investigate  the  first  transverse  mode  for  which  SnU  *  1.64129,  and 
then  Iho  values  listed  In  Table  III  follow  directly.  Note  that  tor  all  these 
cases,  h,  ,  Is  Identically  zero.  Those  results  are  plotted  In  I  Igures  14.2 
and  14.3, 

If  wo  examine  the  curve  of  %  versus  (aJ  for  chamber  number  I 
(  «  0.1),  wo  500  that'  this  curve  does  not  exhibit  a  minimum;  however, 

If  the  values  for  ViL  and  are  extrapolated^  out  to  a  value  of  (a)  «  1.50, 

then  there  Is  obtained: 


U) 

hi  Re. 

hi 

hjRa 

% 

£ 

1 .60  ” 

0.'  148 

0.0172 

-  0,55 

"  0 . 30*** 

3.00 

2.97 

1 . 50 

0.155 

0.170 

-  0.46 

0.435 

3. 14 

3.46 

*  A  qualitatively  correct,  but  not  rigorous  procedure. 
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and  il  Is  ofeservod  Mial  ft-  takes  on  its  minimum  valve  of  approximate I y  2.99 
in  the  vicinity  of  44*  1, 65. 

Examination  of  f  Iguro  14.2  also  shown  Dial  I  ho  curve  of  'Yb  for 
chamber  number  2  exhibits  a  minimum  of  2.44  in  the  vicinity  of  iu  •*  1.84. 
Ilonco,  a  comparison  of  I  ho  results  obtained  for  'these  I  wo  chambers  scorns  to 
Indicate  that  when  the  Chamber  exhaust  velocity  is  larger,  the  stability 
limit  is  small  or,  which  moans  that  the  susceptibility  of  the  chamber  to 
combustion  instability  is  then  greater.  (See  I  Iguro  14.3)  Wo  further  not© 
that  the  minimum  f l  occurs  at  a  different  frequency  for  each  of  those 
chambers*  That  is,  the  chamber  with  the  higher  exhaust  velocity  (at  the 
nozzle  entrance)  lias  a  higher  neutral  frequency  for  minimum  9t  > 

A  very  Interesting  result  Is  obtained  in  an  Investigation  o( 
transverse  mode  Instability  when  duo  to  some  peculiar  combination  of  events 
in  the  chamber,  the  complex  quantity  V)  ^  vanishes  identically  at  a  pi  von 
chamber  frequency.  For  such  n  situation: 


2,-/1.  (l  -  6,’^ 


so  that  either 


n  ~  o 


uRe  (  u  C0'SuJ  ^  ~  ^  1 in  u.2  f  ~o 

V\  i  (i  *  co$ui  £ )  +  ^  l  Re  sin  w  h  1 1 


or  both  hold  true  sjmul'l'aneousiy.  Such  u  solution  involvns  a  certain  amount 
of  indeterminacy,  however,  some  general  conclusions  may  nevertheless  be  drown. 


-  ■  '  u 


is  i  ndetorm  I  no  to.  I'ho 


If  II  qua  I  ion  14,7  holds  alone,  I  linn  > 
stability  IliriH  Is  then  given  by  I  ho  hoi  l/onlal  lino  In  the  'K  ,  & 

piano.  Since  I'ho  ''7b  of  any  given  propel  lanl  combination  Is  gonora  I  I  y  finite, 
fills  would  Indicate  Hitrl  combustion  Instability  can  exlsl  al  that  frequency 
for  which  hj  vanishes,  regardless  of  I  ho  value  of  tho  sons  1 1'  I  vo  1  lino  lag. 

Now,  Ihe  sinu.il  I onoous  solution  of  Eqs.  14.8  ylnlds: 

M.9 

(K  «  0,  I,  2,  3,..,) 

and  In  this  Instance,  'Yb  Is  indeterminate.  Tho  stability  Mini  I  may  then  bo 
Interpreted  as  the  vorllcal  lino  &  ~  O  ,  In  the  %  ,  \  plane.  This 
solution  Implies  Thai  If  a  given  propellant  has  a  zero  sensitive  time  log,  then 
combustion  Instability  will  exist  at  the  frequency  at  which  nj  vanishes, 
i  ogardloss  of  I  ho  value  of  the  Interaction  Index  71  .  Thus  wo  conclude  that 
the  conditions  that  make  tend  to  zero,  promote  Instability. 

If  the  behavior  of  In ?  and  Vli,  >  of  frequencies  above  or1  below  ihe 
UJ  for  which  vanishes,  Is  such  that  71  and  &  are  both  negative,  then 
the  chamber  will  be  unstable  for  those  frequencies.  This  follows  from  the 
dot  In  It  Ion  of  tho  stability  boundary,  since  when  tho  computed  o I gon values  aro 
both  physically  unobtainable,  tho  amp  1 1  f  lea  I*  Ion  cool  flc  lent  A,  Is  then  flnlto. 
I!5.  Compar  I  son  w  I  th  ffxper  I  merit  and  Cone  I  us  jpns 

In  fills  section  we  shall  endeavor  to  compare  tho  theoretical  result's 
presented  In  tho  main  body  of  this  thesis  with  pub  1 1  shod  expor  I mental  data. 

Wo  have  developed  a  theory  for  a  complex  physico-chemical  phenomenon  hosed  on 
a  hypothetical  model  of  the  combustion  process.  Tho  Justification  for  such 
an  approach  resides  In  the  fuel  that  nol  too  much  Is  known  about  the  myriad 
factors  which  Influence  tho  behavior  of  the  entire  system.  However,  tho 
validity  of  our  approach  can  bo  assessed  only  In  terms  of  a  d  I  reel  comparison 
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between  I  hi?  trends  prod  I  c  tejcf  by  1  ho  Theory  and  t  ho  results  obi  o  I  nob  I  o  f  i  om 
expor  I  moot  . 

[ho  pi  tlnary  cons  fdorat  Ion  whltli  proven  Is  such  a  CM  reef  compos  Ison 
Is  the*  foci  I  ho  I  111©  thoorol Ico I  Ireolwnl  presented  here  is  for  I  ho  stab  i  II 1 1  y 
of  a  rock el  system  against  d I s I urbanuos  of  small  amplitude  ( app I  I  cab  I o  lo  the 
on  Sot  of  combustion  I  ns  tub  1 1  it  y) ,  while?  the  majority  ol  pub  I  I  shod  results 
deal  with  I  ho  fully  developed  Instability  character  I zed  by  non- 1  Inoar  effects 
including  shock  waves  and  viscous  damping. 

Because  this  situation  was  recognized  by  Croce©  some  tlrno  ago,  an 
oxpor imontul  program,  supported  by  the  Bureau  of  Aeronautics,  Department  of 
the  Navy,  has  boon  under  way  at  tho  damns  Forms  la  I  Research  Center,  Pi  Incelon 
University,  for  the  specific  purpose  of  studying  the  development  of  combustion 
Instability  In  liquid  propellant  motors.  A  comprehensive  treatment  of  tho 
experimental  findings  is  given  In  tho  recent  work  of  Matthews  (Ref.  23}* 

Since  I  ho  major  part  of  Matthews1  work  deals  with  the  experimental  dot  or - 
mi  nation  of  I  ho  combustion  time  lag  In  motors  operating  at  essentially 
constant  mixture  ratio,  with  a  modulated  propellant  Injection  system  which 
produces  ta  low  frequency  oscillation  In  Injection  velocity  and  propellant 
flow  rate,  we  cannot  make  any  direct  comparison  with  his  work.  Additional 
experimental  effort-  along  the  lines  of  tho  theory  developed  hero  is  clear  ly 
inquired,  before  such  a  direct  comparison  Is  feasible. 

In  any  case,  we  must  look  elsewhere  for  a  verification  of  I  ho 
analysis.  The  accurate  determination  of  quantities  oscillating  at  high 
frequencies  Is  always  a  difficult  undertaking,  however,  aside  from  flow 
visualization,  the  quantify  which  is  perhaps  most  easily  determined  Is  the 
frequency  of  chamber  pressure  oscillations.  Thus,  our  theoretical  results 
may  be  compared  with  experiment  with  regard  to  two  separate  measurements  or 
observations,  rho  gos  motion,  and  the  chamber  pressure  frequency  during 
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uns feady  molar  operation. 

I  pi  us  begin  with  I  ho  mo  i  Ion  of  a  gas  par  Mole.  II  Is  of  some 
Inlerosl  la  first  sketch  the  modes  of  oscillation  at  an  axial  station  for 
the  values  of  5v»U  given  in  f’q.  5.213.  This  may  bo  accomplished  by  noting 
that  al  an  axial  station,  the  perturbations  In  ges  velocity  and  pressure 
may  bo  written: 


Ve.'  ~ 
Yr 

Ye'  ~ 

p ' 


3n  (Snhr)  Coin 9  e,lAJ 

djjl  (SnV^r)  co$n9  e'UJ 

dr 

3n  (sninr )  n  VT0  £ 


r- 

Jn(^hr)  C 


The  sketches  shown  In  rig.  15.1  then  represent  the  Isobars  and  the  Instant¬ 
aneous  directions  of  the  gas  particle  motion  In  a  tangent  la  I -rad  la  I  plane, 
at  a  given  station  %  at  time  *t  .  In  order  to  consider  what  transpires 
as  the  timer  Increases,  we  observe  that  the  resul  t  depends  on  whether  we  have 
standing  waves  or  traveling  waves  (or  both)  In  the  chamber.  Standing  waves 
may  be  I  dent  I  Mod  by  I  ho  presence  of  stationary  nodes  or  nodal  diameters, 
Indicating  that  two  trains  of  waves  of  equal  amplitude  and  frequency,  but 
out  of  phase  by  180°,  have  traveled  past  each  other  continuously.  Such  a 
situation  Is  possible  only  by  I  he  process  of  wave  reflection  at  the  boundaries 
of  the  rocket  chamber. 

The  standing  wave  form  Is  easily  obtained  by  noting  that  Increasing 
time  progressively  reverses  the  pattern  of  particle  motion  and  changes  the 
algebraic  sign  of  the  excess  pressure,  so  that  when  cut  ®ln  ,  I  ho  motion 
is  completely  reversed  and  a  pressure  deficiency  exists  where  there  was  an 
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1 

i 

•a 


excess,  and  vic»-vei  sa.  At  U/\  2Jfl  Y\  ,  tho  ut  I  y  I  rial  pfelufe  is  restored. 

In  1 1 1  o  Iraveling  wave  form,  or  spinning  form,  I  ho  pal  terns  shown 
in  I  ig.  1  *>- 1  con  rotate  continuously,,  since  l-qs.  I  b.  I  con  olso  be  wrillon  in 
I  ho  form: 


Vi'  ,v 
\/  '  • 


Vo 


\ 


'V 


p 


A./ 


,  fwi  i  n0) 

Tn (dhr )t 
dr 

JjAi±r)  g 
r 

— r  /  i  \  t’fwi  +  nt?) 


I  ho  1 1  mo  for  a  comp  I o I e  r ovo | ut I  on  1 s  g I  von  by 

II  we  restrict  i  r  attention  to  the  I  I  ml  tangential 
mode,  and  nolo  I'hot  I  ho  particle  velocity  Is  superposed  on  a  moan  gas  motion 

V-  ,  the  resultant  par  I  I c; I o  mol  Ion  may  bo  sketched  as  In  Fig.  1 5.2,  In 
which  the  gas  propagates  axially  as  It  spins,  so  that  the  motion  of  an 
Individual  gas  particle  Is  somewhal  Ilk©  a  corkscrew.  The  numhor  of  Complete 
devolutions  that  an  Individual  gas  particle  makes  In  I  ho  chamber  can  be  deter- 
mi  nod  exactly  by  considering  the  Lagranglari  derivative  of  the  par  tide  mol  ion, 
howovoi  ,  an  approximate  value  Is  given  by: 


and  I  t  Is  clear  that 

At .  . 

w 


N 


WJ 


L 


wn  V*  avk„ 


Note  that  if  there  were  no  mean  motion,  as  In  a  cylinder  with  closed  ends, 
all  the  gas  in  the  chamber  would  either  "slosh"  or  "spin"  slmu  Maneously  at 
every  axial  station.  It  Is  also  observed  that  the  linearised  analysis  ponnils 
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the  simultaneous  oxlslenco  of  any  combination  of  standing  and  traveling 
waves  w i II i  arbitrary  direction  of  relation. 

In  Iho  expor  Imon  la  I  I  rives  I  I  <ja  I  Ion  i  o  ported  In  Ref.  I  b,  probe 
microphonos  wero  ulllUod  lo  dotormlno  Iho  frequency  and  Iho  amplitude  and 
phase  of  Iho  pressure  oscillations.  On  the  basis  of  llioli  measurements,  Iho 
authors  claimed  that  the  first  Irons verso  mode  could  ho  do  toe  tod  In  the 
standing  wave  (sloshing)  form.  On  the  other  hand,  streak  photographs  hove 
boon  lekon  through  Iransparonl  sill  windows  by  other  Invest  I  gators  (see 
for  example  Rof.  21),  which  Indicated  Iho  prosonco  of  a  rotating  luminous 
zone  which  propagates  on  a  helical  path  along  the  chamber  length.  In  this 
particular  case,  Iho  authors  state  that  Ihelr  experimentally  determined 
frequency  Is  approximately  that  of  the  first  transverse  mode. 

those  results  give  credence  to  the  possibility  of  obtaining  bath 
standing  and  traveling  wave  forms  during  unsteady  operation  of  the  com¬ 
bi, is  I  ion  chamber.  Our  analysis  anticipates  the  results,  obtained  by  froguoncy 
measurements  since  wo  have  already  observed  that  the  frequency  of  neutral 
oscillations  UJ  as  determined  by  an  actual  solution  for  Iho  eigenvalues 

lA)  and  £,  ,  will  have  the  same  order  of  magnitude  as  the  value  ol  uJ 
given  In  hq.  7.6  (Iho  exact  solution  of  Iho  wave  equation  In  a  cylindrical 
chamber  with  closed  ends).  I  t  Is  fm  flier  noted  that  since  the  experimental 
value  ol  the  frequency  dur  Incj  unsteady  operation  is  very  nearly  given  by 
I  q.  7.6,  some  Investigators  have  assumed  that  Iho  gas  dynamical  behavior 
of  the  chamber  can  he  adequately  described  by  the  classical  wave  equal  Ion. 

This  assumption  obvlr  I y  oversimplifies  the  actual  state  of  affairs. 

In  conclusion,  v e  remark  that  since  the  character  1st Ic  equation 
of  a  general  rocket  system  has  now  boon  solved  Implicitly,  for  I  wo  Important 
mechanisms  capable  of  producing  linear  combustion  instability  In  the  high 
and  Intermediate  frequency  ranges,  this  solution  may  now  bo  used  as  a  research 


f@ol  In  wlii<h  a  Dr  oafl  i  n  vos  I  I  ga  I'  ion  of  iho  behavior  of  different  roe kef 
mold  s  is  niacin.  In  such  a  paramoti  ic  Investigation,  I  ho  susceptibility  of 
ilia  rock  of  mold  10  transverse  wave  or  onlropy  wave  In  stab  I  I  I  ly  can  be 
determined  for  d  I  f  for  on  b 

I)  dlr.lr  Ibui  lono  of  combustion 
/)  chamber  goomoliy 
■  3)  I  ejection  systems 
4$  t-xho.usl  nozzles, 

In  view  of  the  comp  I  lea  I od  form  of  Iho  solution,  I  I  Ik  unlikely 
that  the  general  behavior  of  a  Systran  will  bo  os  I  ah  I  I  shod  without  a  major 
of  fort  along  the  linos  of  additional  numoricul  Computations.  Iho  Iwo 
numerical  cases  floated  here  ropro&onl  two  ol  Iho  morn  Interesting  results 
which  liavn  already  boon  obtained  In  III  is  program  of  investigation.  ©real 
difficulties  worn  encountered  during  the  course  ol  performing  the  computation 
Negative  resull'S  (un  intorosl I ng  values  of  the  eigenvalues)  were  obtattled  in 
a  number  of  cases  which  war®  therefor©  not  Included  in  this  presentation-. 

It  Is  fell  that  I  hi  t  wo  cases  wh  I  ch  are  I  ric  I  ruled  I  ft  III  I  s  ana  I  ys  I  s  are 
typical  of  the  results  which  may  be  obtained  In  an  ana  ly  I  leal  In  vosl  I  gat  Ion 
of  combustion  Instability. 
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Suporgr  !  I  I co I  Oscillatory  Dlscharqo  of  Convorq I n < | ~D  I  vorg i n p  No^ ^  1 05 


Ilia  oscillatory  discharge  of  convorgjfiq-d iverg I ng  noz/lo  ope rat I ng 
in  i  ho  *j  u  p  o  r  o  f  it  tea  I  range  has  boon  treated  by  Tr.  Ion  (Hof.  6)  ant!  Crocco  diets, 
olid  19). 


Tslon's  treatment  was  one~dlm@nslonal  and  was  rosti  I cted  to  the 
case  In  which  the  oscillations  In  tho  Incoming  I  low  were  i  sot  hoi  inn  I .  Further¬ 
more,  his  solutions  wore  obtained  for  very  low  frequency  and  for  tho  asymptotic 
case  of  very  high  frequency.  Of  those,  the  latter  of  the  two  Is  beyond  the 
range  of  applicability  to  rocket  motor  oscillations^  while  the  former  Is  not 
sufficiently  general,  Therefore,  Crocco  ox  ten  tied  this  treatment  lo  tho  general 
coso  of  longitudinal  arid  transverse  non-l sothenrw I  oscillations,  and  he 
obtained  $r 'utlona  over  the  full  frequency  range  of  Interest* 

We  have  used  tie  results  obtained  by  Crocco  ns  tho  houndary 
condition  for  the  combustion  chamber  analysis,  and  have  wriiton  these  In  iwo 
different  forms.  For  tho  analysis  of  transverse  per+urbat Ions,  we  have  taken; 
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7.7 


whoro  and  C  are  the  nozzle  admittance  coef f lelon+s  and  are  complex 

functions  of  tho  frequency  lA)  ,  Swlo  ,  and  the  nozzle  geometry.  This  equation 
states  that  for  a  given  nozzle,  the  perturbations  at  tho  nozzle  entrance  must 
bear  a  certain  relationship  to  each  other  with  regard  to  phase  and  amplitude 
If  neutral  oscillations  are  to  be  maintained  for  a  given  mode  and  frequency. 

For  tho  purely  longitudinal  case,  wo  have  taken; 
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II  Is  clear  that  when  Snh  °  0,  we  are  treating  the  purely  longitudinal  case 


•'rid  lq.  7,7  mu  si  reduce  toEq,  12.14,  for  ono-i  imore*  Iona  I  flow,  t  ho 


subscript  ?•  i'  no  longer  necessary  and  them  wo  obtain; 

./I  -  -  oTv, V«  5;  C  -  -^Ve  a  i 

II  Is  of  some  Interest  to  con-  I  dor  the  admittances  tfvi  and 
In  greater  detail  In  Reference  6,  Tslon  has  shown  that  for  nozzles  which 
have  o  I  i nodi  variation  of  velocity  In  tholi  subsonic  part,  1  ho  governing 
equations  Is  a  non -homogeneous  hyporqoomotr  ic  differential  aquation  In  the 
form: 
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go  that  at  the  no/z lo  entrance 
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and 
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non -d  I  mon  9  Iona  I  frequency  ^  which  may  be  related  to  the  uj  of  th« 


chamber  which  appears  In  all  the  chamber  equations,  and  hence  UJ(^  will  ho 
used  as  the  Independent  variable  for  all  the  curves  of  the  admittances. 

The  relevant  condition  to  ho  satisfied  of  the  nozzle  entrance  is 
given  by  Crocco  (Rol.  7)  as* 


cPe  ~  X  0*  e  ^  A  S 

In  this  same  reference,  a  series  solution  was  obtained  for  low  frequency  by 
taking:  ,  .  .  .  \ 

rfo)(^)  +*/2cr  1  ($ )  0  f^yi . 

B  Zio)  (') )  +  i(iz>l'H$)  +■  (</r  . 

however,  the  series  was  cut  off  after  only  two  forms  and  the  results  wore 


correct  to  0(  c^3  ).  Since  that  reference  first  appeared,  some  exact  numerical 
solutions  have  been  obtained  over  a  wide  frequency  range,  and  direct  comparison 
shows  that  {iv\  and  fty\  may  be  determined  analytically  with  sufficient 


accuracy  «jn&n  iaJJ  <0,3.  by  lie*)!  cm  !  cu  la  ii  ncj  an  add  i  Maria  I  term  In  the 

3 

sorie'..  Wlir>n  Ibis  ir,  done,  there'  resiilts^ 
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For  example,  taking  K  «  1.20  and  Vt  0.10  ,  wo  obtain) 

iYm  -  ,|  [  |  4.  ZZ,lOOt//2  +  6  ^ 

fhn  *  .&[  I  -  Z,  ISO  1/3  -  j 

As  stated  previously,  all  graphs  of  t(\r\  and  f2/\  are  plottod 
utilizing  as  the  Independent  variable.  This  quantity  Is  related  to  t ho 

frequency  in  the  chamber  by  means  of: 

urfa  _  &  _  rV  "T+T' 

(aj  "  C>  k  L*  ^  ^ 

where  K  Is  I  ho  non-dimensional  velocity  gradient  In  the  subsonic  por  tion  of 
the  nozzle: 
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Thu s ,  I  ho  smaller  Ik  Is,  the  longer  Is  the  subsonic  portion  of  the  nozzle, 
Since  the  length  of  tho  subsonic  portion  of  the  nozzle  fs  given  by: 
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Thu  graphs  of  (ft  ,  w  mid  Caro  likewise  plotted  versus 
In  this  cose,  the  frequency  In  tho  chamber  (a)  Is  related  to 

[Ki u  „ 

^  ~  Tr? 
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by: 
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since  another  scheme  of  non -d I  mens  Iona  I  I zat loo  was  utilized  In  these 


computations. 
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APPENDIX  U 

Dor i vnt Ion  of  Injuclor  Response 

Wo  will  determine  hero  the  variation  In  mixing  ratio  and  Injection 
velocity  as  a  function  of  the  or.cl  Mating  prossuro  uj  fho  I h J oc lor  face. 

Th  I  will  bo  dono  by  flrsl  cone  I  during  I  ho  frequency  response  of  a  single 
food  I  I  no  and  t  hon  so  looting  tho  don  I  rod  Impingement  pattern  to  yield  I  ho 
combined  effect  of  propel  Icirvl  and  oxidizer.  The  rolevairl  equations  for 
quas  I  -orio-d  I  men  s  I  oria  I  flow  In  I  ho  food  I  I  no  are: 


a  ^ 

which  are  the  conservation  of  mass  and  momentum  respectively,  and  where  A 
Is  I  ho  cross -sect  Ion  of  the  food  1 1  no.  Tho  coordinate  system  Is  Indicated  In 
fig.  m. 


Upon  non-d 1  mens  Iona  1 1  zing  as  follows: 


Introducing  small  perturbations; 

V,  (et)  =  Vx(b)  +-Vjt .'(*,*) 
p  C1!*)  =  }'  (■■*)  -I-  ]>'  (*;t) 


con  i  I  nu  I  I'y  and  rnomon  I  um  y  i  o  I  d: 
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Wo  may  separate  the  variables  by  taking: 
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B  Ic 

u  y.b 

D  2c 


V/M)  *  1(^)e 


B  4 


arid  then  wo  obfaln: 


o 


s* 


1 


,  if 

•A  d* 


B  id 


e  2d 


Hof on  solving  these  aquations,  I  of  us  briefly  consider  the  boundary  condition 
A I  the  Injector  loco,  ( %  0),  tho  steady  state  pressure  and  Injection 

velocity  have  tho  I  r  <’  sign  values  while  the  pressure  perturbation  must  be  the 
eluant)  or  pressure  perturbation,  and  tho  as  yet  unknown  velocity  perturbation 
will  correspond  to  an  Initial  condition  for  the  chamber.  Because  of  tho 
re  I  at i vo I y  small  d I  a motor  of  conventional  injection  ports  with  regard  to  the 
reserve  I r  proper ,  wh  I  ch  means  that  on  I y  a  sma II  liquid  I nert la  Is  I nvo I vod , 
and  because  of  tho  further  likelihood  of  entrapped  vapor  In  the  reservoir,  It 
seams  reasonable  to  assume  that  tho  reservoir  pressure  Is  constant  only 
a  short  distance  from  tho  port.  That  Is,  the  reservoir  proper,  Is  an  Infinite 
source  of  liquid  and  the  velocity  vanishes  whore  the  pressure  Is  constant. 
Hence,  we  may  summarize  tho  boundary  conditions  ns: 


a*o;  Ve  (°)  = 

-  to 

2=Er-  YtM=° 

g  (2-p.)  =  o 


pM  =r 
W°)  --  % 

pftp.) 1  Pit. 

'P  (  4  ft )  =  O 


B  5 


integrating  £q,  I.1  2b  and  Applying  l ho  boundary  conditions,  wo  obtain: 


o  - 


+  z 


and  B  lc  yield: 


fa  h 
fa  &  ^ 


YY  i  ~  Con 


fY\  =  co  n 


=  const.  =  _JiL 

yJO  vio 

Now  Integrating  Eq .  13  2d  from  to  *t5  «  0,  wo  obtain: 


n(*)d 


2:  4 


x  n* 


On  Introducing  Lq.  D  /,  I  ho  first  term  yields: 

li  1<'|Jy  V*  C  m 

Uj|  W.  r'/’il  I  rr.  F*  (  °deT  I  n  Y\  i  i  1  uPi  1  f’i  41  I  k 


p. 


where  ^  a 


a  an  equivalent  length. 


Substituting  back  and  solving  for  v]c  ,  we  find: 


whore  the  quant ity 


xfa  YeoO  +  ^) 


a  = 


Vjto  Ldf. 

a  dimension  loss  characteristic  Nino  of  the  food  i  I  no. 
Wo  ma y  d  o fine  I h o  I r  a n  s f o r  f u n  cl  1 o  n : 


*  TO  V/i  %o  0+$a) 

which,  for  neutral  oscillations,  becomes: 

_  »  b  4"c\  vo '  eu  <\, 

q(w)  .  - - 

Xfi  VjlA  )+  T’lX1- 


xfa  Viol  \  +  u)^a v- 

Arid  now,  let  us  use  this  result  In  der  iving  the  response  of  a 


typical  hi  propel  I  ant  Injector  where  we  will  let  the  subscripts  ox.  and  f 
denote  oxidizer  and  fuel  respectively.  If  we  consider  an  Injector  of  the 


typo  shown  In  Fig.  \\2,  then  tho  longitudinal  component  of  Injection  velocity 


V/,o  ■=  (tyo^yO.x  +  (Vio  CoS^)t 


whoro  wo  liavo  assumed  conservation  of  the  axial  component  or  momentum  during 
Impingement  in  the  steady  stale.  And  for  unsteady  flow,  under  tho  same 


ossumpt ion, 


(uj)  Co 4'  1$  )oX  4-  o)  Coi 


whoro  we  have  assumed  that  tho  fuel  and  oxidizer  ports  experience  the  same 
Instantaneous  chambor  pressure. 

Proceeding,  wo  lot  F*  denote  tho  mixture  ratio,  then 

r  =  P'  +  r1  •  +  B  16 

VY'-P  rn  -f  (ity^-f)  ^  ^ 

where  wo  have  noyloetod  second  ardor  terms  and  where  na  J?2r  denotes 

/ 

tho  fractional  perturbation  In  mass  flow.  Now, 


Mi)*  ^  =  HEEll  =  Ei.  - 
Wo  ftv^Ao  Vi0 

end  hence  substituting  Into  Eq.  B  16  we  obtain: 


if  p 

Via 


Vg„  /ox  l  V/o/f 


B  in 


Wo  may  now  obtain  a  relation  for  the  fractional  perturbation  In 


tho  Injection  rate.  Since, 


me  -  me  4  -  mox(  14-^oxJ  4-  ) 


there  fol lows 


VYU  _ 


ww* 

I  4- 


JK  4”  - ss 

I  +  p*  ^ 


and  on  Introducing  Eq.  B  17 


•  1  I 

<>nd  wo  nolo  that  although  B!  111.  for  AnHl  .  ... 

'0f  Gac*'  t  ood  1 1  no ^  wo  cannot  put 

">•  B  I'.  Into  the  some  form  ns  Eq.  0  20,  since  In  general  (co t fi)  ,  ond 

(cvs/i)$  r re  not  equal.  Thus 

/Wt\  _  f  ($Mc°.s/3)o*  +  (aMcoj/S)|  h  .i 

WJ  l  (Vi.co^).,  +  (vfc.COJ/S~  j^e  1,21 

Honco  the  transfer  function  for  the  case  of  »  blpropo, lent  Injector  Is  given 


whom 


jM  = 


■jf  =  J  M  vxo 

=  f  CoSA),k  +  (  ywJtojA 

\  —W»IIIM'./<hwii r~lwiu«  I  |wtIM|^|()[ i,tkT  l 

(  (yiucoj ^3  )oX  +.  (V^cos^)^ 

We  may  summarize  those  rooults  by  writing: 


r  )„  ■  GH%ed 
jL)o  .  H(wjv?e5t 

III  = 


“l,Sr0  Jnd  aro  tloflno(i  through  Eqc,  B  10  and  B  20.  Wo  note 

that  If  another  injection  system  Is  utilized,  then  the  analytical  form  of 
the  In Joel  or  response  Is  nevertheless  given  by  equations  of  the  form  B  24. 

One  additional  quantity  will  bo  derived  hero.  *  have  stated  else- 

I  I 


where  that  wo  will  assume  that  „  0,  0„d  wo  may  therefore  obsoryo  that 

each  particle  will  retain  Its  Initial  value  of  h/,  ,  kh  being  a  constant. 
This  enables  us  to  wri to: 

WsM  =  In/,  (o,  i-(*-±L-  \ 

v  '  )a  Vt  li'Mv)]  I 

.  ,1  ,  ,,  B  25 

aV,tr.  t-/‘ 


V/(a')  > 


*  J  ?J 

°  Yite1) 


8  “6 


I  ho  spec  iih  enthalpy  Vt/n  of  the  propel!  ants  depends  on  the 
Instantaneous  value  of  the  mixture  ratio  and  hence  the  enthalpy  will  ho  larger 
or  smaller  depending  on  the  mixture  ratio  which  In  turn  depends  on  the 

pr  v‘Uf0  Qt  +h®  cht,ml,or  face.  Expanding  \y  In  a  Taylor  series  about  Its 
steady  state  valuo,  wo  obtain; 


and  doflnlng 


h  I  *t 


ihi  (k 
ar  v 


r)  +• 


iK 


ht« 


( ihi 

\  d  K  /o 


we  obtain 


iK 


Thus  Eg.  0  23  becomes; 


=  aK  K(  )  ft"  f-V/T'll  +  (»-') 5  *7  Wo1  \  ri  /tj*1 

^r/0  WJ.  -*•  ( 7X  J  Lt  1°  %(» 


r  /o  "■  w 

and  substituting  from  Gq,  n  24  wo  obtain; 


B  26 


B  27 


0  \(v) 


,1 


B  20 


*  M(w>) 


B  29 


wh  1  ch  dot I  nos  M  (W) . 


C-l 


Append  Ik  C 

Dor  I  vat  Ion  of  Equation  13. 

.  rr-Ti ' ....  m * rnf  r 


W®  Wl11  bogln  bV  Introducing  Eq.  12.15  Into  Eq.  12  16,  and  them 
•oltlng  X  ( it '.)  -  <f  (constant ) ,  thoro  follows: 

( ^)f>  r  (i  ~  e~tw5 

+  ^  i  H  |  W  )  -  £  t  UJ 


cosco  z  1  dy  (z1) dv 

ctZ' 


dr 


■>V{WJw-i5inw  +  lMJ'[vEc%<)cotM(i.ii)'jjii  c 
[  If  F(0,(2')  Smwf|_  gi)  Jt/?1 


f'  I 

/o 


which  may  bo  rewritten  as 

(•) 


o 


(  5  IK1  (aJ  —  ^  i4  cod  UJ  Jr 


-  i 


+  lUc^M  +■  B(w) 

whoro  wo  havo  Introduced  the  notation 

If  ~  V  ^H(w)-7H-T(wj(; 

LM 

1  W 

v  yt  C :  ~  e  -‘wf  j 

and  upon  Integration  by  parts,  It  may  bo  shown  that 

D(w 


pMUtMSMI/j 

H 

I — ..I 


(  /  if  (i1)  Jv  (gi)  cj 


coiwh-a'J  coi  wi1  dy 


Ve  cojcu 


/«£> 


si?  (*')  HO  /i  -  2tL()  J 
o  oil*  1 


C  2 


C  A 


1,1 1’  '  Ma/)  consists  of  the  terms: 

Bfw)  =  u»  |  (o-t)  j  v/i1) coiiwe1  sinw 


n,n]h<) 


*  Sintu(l~2’)cl2 


UJ  Jo 


!51Wi  j 


i')  sin  we  1  5  in  wp  -  2()d% 

a  I  0  1 1/1 1.  I  /  . 


-ft'Jrfa 


f  (a-y)  jj  V»  cofw(|.e 

+  (3-(j')  jo  V^1)  Sinwja  1  coicu  (i~2i)c/ 


&  ')c/r 


4-  l  e^nM 

Ml'  >i  >  foul  th  and  seventh  terms  may  be  expanded  and  then  combined  to 

y  w  id*  l  '  | 

4  ”  wSlflwJ)0  y (*0ti * 1  -  (i-v)w  j#  v(t'\  sin wfri Zi'jda 

(L  f  I  j  i  ^ 

jT  )o  fa  'b'  r'°  s  U)  ■  CoS  oU  (  I  "  2  B1)  |  c/ ?. ' 

It  wo  now  define i 

/Un  ■  -*U)  #W  i" p 

A  r/  s  .i 


ihen  from  £q.  12.12  wo  obtain: 


Vt  (>'j 


^0>4)  =  -JHe^^Afa) 

Ju  d  ( ro^AvN  -T,  v 

%  3i  $()  *  Ji,^c  #(*)- 

V  to  /  di 


/If*)  Cut  j 
JMJ  C  9 


s  i  nee 


4  =  a  fi  a: 

'  '■  d-Z  )o  v^fe’) 


so  that  upon  Introducing  the  notation 

B,(w)  =  /J  foiuj(l-22')d4' 

B,  fw)  -  |o'  Snow  (i-i')da'  I*  e"‘wT^«")  d£i|('a")  c/a" 

Bj  (w)  -  I'  VM  sinw(i -zs'jola1 

_  C  10 

P4(w)  -  f'  i¥(i\)  coiw(i-2')d?,' 

0  d^1 

B5(uj)-  |P  /H2-')  sin  vaj  fi- «')dz‘ 

B6(0J)  -  /J  A(l')  M  S«w(|-B')di' 

B7  (w)  -  |o'  AM;e’lwT*^  V(s')  coauj(i-.a')d«l 
Bj(w)  ■  JP  M.  (a’)e'",w't(^4^tJinujp-a')dal 

—  _  Q  £ ' 


Eg.  C  6  becomes; 

£  s  ujsiniA)  f  -(z.-Jr)uj  E>s(w) 

T  B  /Q 

jU 

z 


and  then  Fig.  C  5  may  ho  written; 


B Ml  =  -  ifooBilw)  H  ( vaj )  4  VuABjw)  G(w)  iwf 

4.  J(w)  ~  2u/Bs(ia)^  +  ixfe.’E»4(w)  +  c  12 

-(L&»HJ  +  5fUj /M|w)  ^4(uj)  ' 


*'*'’4«wenuk*.»S. 


■'-’-f  „  rH ■*■*»  *********  ff )  tffitu "•**■**  , 


04 

W°  "My  1,01,1  ht|'  12,17  "'  u  sln,ll°'  mornior.  first  substitute 
'  ’■■>>  Into  Eq.  12717  so+  ?&).  *  (eonstant)/  Tho„  wo  obt#|(|  ^ 

result: 

/✓i 

coj)  uj  4*  2  l  sin  uj 

1 


n 


+  Mo  |  irE^^i')  sinuj(i-gi)  !  dgi 

-  IW  |o  jjf  |  M1)  coj tAj^ I- e«)  I  ^ 

which  may  ha  rowr  Itton  as; 

/  ^fc 

(  ~0Q  )  ““  ^  4“  2  c  Ve  5mu.j 

~  1  2  C(w)  4  L  A(w) 

whoro  upon  Integration  by  porta,  It  may  bo  shown  that: 

CM  =  wJJ  vfH1)  eojujf,.^)^ 

Lml  AM)  consists  of  the  terms: 

•,  ~~  "*  '  -  "ft"  ‘  Kt 

AM)  ~  uj  V(%!) "coiuj fc'cosw  (1- a ,1) 

-  IT  ['  coswfi-e'Jdt'f^'e-^^M  dy 

I  '  (q)  '  d*" 

*  t* 

*  ^  }0  fx  (*')  S,nuJ£'  CoSujfl-fci)^! 

~j  (o'  costv(i-e')  da1 

(3  C  7ft')  ^E')  sm  -u^i-  g'jcM1 

|J  V(t') ..sintw-a’  sin o>(i~ 2  ')da' 

-  *M(w)  |J  n^Ttfa') 

(.■  2 1 


c 


■«  1  /. 


c 


C  lb 


C  16 


rlw«.  -^JKiSh; 


I 


1 


The  first,  fourth  and  seventh  forms  may  bo  oxpanded  and  then  coniMnod  to 


yield* 


-fA  =  wCoiw  jj  -  (l-H )wjj  W*')  COSUi(,-Z«')Hi' 


•t-  -4  P  [  snow  -  iiKiw()-2a')] 'd* 


And,  upon  Introducing  the  notation: 


AiM  =  I  '  h(i')SWM  (t-2S')d4l 


Aj.M  f 

A  3  N  = 


co5w(i -£')ck'  f*  dYfa'jck1 

v  *>  ole' 


A3  M)  r  j  V  (%')  coi'tAJ  ( I  -  2e')  1 

f  dy(£))  5imu(!  -  t£')oI£ 1 

^  dxL 

Ast^)-  ('  A(i')e'lVj^^^  Cojwfl-  *')  d  s 1 

Jo  ^ 

h(*  (W)»  jn'  A(i')e'l'wTt(v)  £(*')  coSw(i-E')<Ja' 

A  7  lw)=  f '  A(*')  v(a')  Sinw  (i-e1)  olz1 

0 

As  (w)j  ('  n  (4.)  coJwfi-e'jcTz1 


0  di' 


Eq...  C  1 7  boeomon ; 


»  UJCO.SW  [J  v(^')d%'  -  (^.-v)uj  Ai(w^ 


+  ^  S.io.to  [o  -  4/A,(>aj) 

while  Eq,  0  IB  may  be  written: 


'  C  19 


AN  ~  -vwAzMH(w)  t  jrw^ N&N 

■ 

*>  f  At 

4  J(uj)  -  2UJ  f\$(w)  +(JlA(o(w)  **  (3-^)  LtO  A7  (to)  C 


-i  Al(w)  ]  -  JfUJ  Mfuj)/\^(w) 


At  this  point,  we  may  Insert  Eqs,  C  2  ,  C  14  and  10.32  Into  Eq.  12.14  and 


wo  obtain: 


\ 


**■*  — "  ■  ■ 


f- 

. 

$ 

m! 

''t-: 


.  '•  •■• " :  •  '•  $  ’W‘a 

M' 

c-6 


. igfiitK 


hiinf  =  k  %(\  -  e'iw^ )  H'  k 


whoro 


k  r  vG(w)[ui1-B2.|w)-wL(w)J  +  isVk  i/nwMtlwjS^J 
^  z  -  -i  Jr  Vt  p<n  C(uj)  -  s  Diw) 


find 


V. 


loi  -  Ve  cx!n  ( C.OJW  +  iAiinw) 


i  i 


c  ?\ 

C  22 


A-Av- 

&\N>; 
w 


'  \ 


' 

wmf 

OfM 


Is*® 


+  Ve^^  )  \  -  ifwAiMH  M  +■  j(w)f 

/ 

Jlo  /  A  .  f.J  )  /»,  v,  Nr, . .  y\  /..,\  'A  /  \ 

-  VUJ 


/k,  {  A  (a  (vJ  )  (<3i- &)iUJ/47  (ui)  ~  I,  ^  ( UJ  j 


“V  Ve/Jn 


li  \y 

ft 


'  \ 

■  if  H(w)  L(im) 
7(uij 


C  23 

(w) 

+  £  SMU)  +  if  V<2  CO^UJ 

B  +  iOAJ  B^fw)  H  (w) 
■2u)6sM  + 

":Vl,/|jJj  J  -  JTUJ  M ftA> )  BafwJ  . 
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